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SOMMAIRE 



Soit D le disque de rayon 1 ou soit S(a) le secteur du disque de rayon 1 et d'angle 
a E (0, 2n). Soit le Laplacien A = + ^ sur A ou A = D ou A = S(a). 
Dans ce memoire, des questions relatives au probleme des valeurs propres de l'ope- 
rateur A sur A c'est-a-dire relative aux solutions de l'equation A(-u) + Am = 
avec des conditions de Dirichlet aux frontieres seront abordees. D'abord nous 
rappelons que l'ensemble des valeurs admises pour A c'est-a-dire l'ensemble des 
solutions pour A sont les zeros des fonctions de Bessel de premier type. Cet en- 
semble de valeurs propres s'appelle le spectre, c'est l'ensemble des valeurs de A 
etant solutions de l'equation precedente. Si le spectre est denote par Sp j4 (— A) et 
si j Vjk denote le k e zero de la fonction de Bessel d'ordre v > de premier type, 
alors 



Le premier probleme aborde consistera a determiner la premiere ligne ou la 
deuxieme fonction propre s'annule dans le cas d'un secteur S(a). L'ensemble 
nodal d'une fonction propre Uj est le sous-ensemble Zj C A sur lequel la fonction 
propre s'annule, c'est l'ensemble des zeros de Uj. Les domaines nodaux de Uj sont 
les composantes connexes du complement de l'ensemble nodal, c'est l'ensemble 
des points du domaines ou Uj ne s'annule pas. II est utile de se rappeler que la 
question precedente ne se pose pas pour la premiere fonction propre puisque cette 
derniere n'a que la frontiere comme seule ligne nodale. Comme il y sera montre, 
si Tangle du secteur est sufnsamment petit, alors la ligne nodale est decrite par 
l'equation en coordonnee polaire r = -%r- < 1 ou v = — sinon elle est decrite par 



SP5(«)("A) 



Sp D (-A) 




; n= 1,2, ... et fc = l,2,...} 
n = 0,1,2,... et A; = 1,2,...}. 
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6 — |. Nous verrons egalement qu'il existe un angle critique a pour lequel la 
ligne nodale est indeterminee et que cet angle critique est donne implicitement 
par l'equation j\{2v) — 32^) = ou v = ^. De meme, lorsque Tangle a = ao, 
nous verrons que la multiplicity de 2 e valeur propre est double. 
Le deuxieme probleme consiste a construire un algorithme pour obtenir les m 
premieres valeurs propres, cet algorithme est egalement efficace pour des valeurs 
larges de m. Cet algorithme ne doit pas resoudre les systemes d'equations diffe- 
rentielles ordinaires relatives au probleme A(w) + Aw = puisque la resolution de 
tels systemes est tres couteuse en temps et en memoire voire presque irrealisable 
pour des valeurs propres tres larges. Le probleme d'ordonner les valeurs propres 
pour lesquelles une caracterisation de deux parametres existe est equivalent a 
connaitre la structure des domaines nodaux des fonctions propres {uj} 1 JL 1 . La 
clef de voute derriere ce probleme est d'utliser un theoreme du a Courant sur le 
nombre de domaines nodaux d'une fonction propre. 

Troisiemement, soit N(X) la fonction de compte des valeurs propres c'est-a-dire 

N(X) = J2 1- 

\j<\ 

Une explication approfondie d'un theoreme du a Kuznetsov et Fedosov montrant 
que 

sera donnee. Pour y arriver, un « grand » detour sur l'approximation uniforme des 
fonctions de Bessel de premier type d'ordre reel strictement positif et d'argument 
positif. 

Finalement, un algorithme permettant d'evaluer exactement et efficacement N(X) 
sans avoir a trouver toutes les valeurs propres Xj < X sera donne. Cet algorithme 
n'utilise que la monotonicite des zeros des fonctions de Bessel et peut trouver 
par exemple iV(9 • 10 8 ) rapidement. Ce 2 e algorithme est nettement plus efficace 
a utiliser pour evaluer N(X) que si nous trouvions naivement toutes les valeurs 
propres Xj < X en utilisant le l er algorithme. 
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SUMMARY 



Let D be the unit disc, S{a) be a sector of the unit disc with angle a G (0, 2ir), 
and let A = + ^2 be the Laplacian. 

In this master's thesis, we will consider several questions arising from the study of 
the eigenvalue problem A(u) + \u = on A with Dirichlet boundary conditions, 
where A = D or S(ct). We proceed by first recalling that the spectrum of the 
negative Laplacian (i.e. the admissible set of eigenvalues A for the above eigenvalue 
problem) is given by the zeros of the Bessel function of the first kind. In other 
words, if we denote by Sp A (— A) the spectrum of the operator —A on A, and by 
j Vi k the k th zero of the Bessel function of the fisrt kind of order v > 0, then 



The first problem consists of the determination of the first nodal line of the second 
eigenfuntion of a sector S(a). The nodal set of an eigenfunction Uj is the subset 
Zj C A over which the eigenfunction vanishes. Nodal domains of an eigenfunction 
Uj are the connected components of the complement of the nodal set i.e. of the 
set at which the eigenfunction does not vanish. Note that this corresponding 
problem for the first eigenfunction is uninteresting because the nodal set would 
consist of the boundary dA itself. We will see that when a is small enough, the 



the second nodal curve is undefined, and this angle is given implicitly by solving 
3i{y) — Ji(2^) = where v — -. Also the second eigenvalue has multiplicity 2 
when the angle of the sector is a . 



SP5(«)("A) 
Sp D (-A) 



; n= 1,2, ... et k = 1,2,...} 
{Jl k ; n = 0,l,2,... et k = l,2,...}. 
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The second problem is to devise an algorithm for obtaining the first m eigenvalues 
that is also efficient for large m. Such an algorithm must not rely on the expli- 
cit integration of the underlying ODEs, since the solution of such ODEs would 
become prohibitively expensive when large eigenvalues are involved. We will see 
that the problem of properly ordering the eigenvalues (for which a two-parameter 
characterization is known) is equivalent to knowing the structure of the nodal sets 
of the first m eigenfunctions {uj} 1 JL 1 . Thus, the key insight behind the algorithm 
is to perform the ordering using nodal set information, which in turn is provided 
by Courant's theorem on the number of nodal domains of a given eigenfunction. 
Thirdly, let N(X) be the counting function of the eigenvalues i.e 

N(X) = 1- 

A,<A 

A detailed explanation of a theorem given by Kuznetsov and Fedosov showing 
that 

A'(A) = ^-^+0(A>) 

will be given. In order to get there, we will look at some results concerning the 
uniform expansion Bessel functions of the first kind of strictly positive real order 
and evaluated at x G R + . 

Finally, we will present an algorithm that allows the exact and efficient evaluation 
of N{\), without having to compute all of the smaller eigenvalues Xj < X. The 
proposed algorithm only relies on the monotonicity of the zeros of the Bessel 
functions and can, for instance, find quickly N(9 • 10 s ). The latter algorithm 
provides a much more efficient alternative to the naive approach of finding, with 
the help of the 1 st algorithm, all the eigenvalues Xj such that Xj < X. 
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INTRODUCTION 



Ce memoire est un exemple de ce qu'il est possible d'etudier en geometrie spec- 
trale. Puisqu'il est generalement difficile de definir et de delimiter exactement 
un domaine de la connaissance, le lecteur trouvera dans ce memoire quelques 
exemples de ce qu'il est possible d'etudier en geometrie spectrale. Neanmoins, 
personne ne se trompe en affirmant que la geometrie spectrale est un domaine 
alliant la geometrie et l'analyse spectrale. Traditionnellement, l'analyse spectrale 
etudie le spectre des operateurs comme les matrices ou encore des operateurs 
differentiels comme l'operateur J^. Generalement les operateurs agissent sur des 
objets comme des vecteurs ou des fonctions. Les domaines sur lesquels ces objets 
sont definis sont plats generalement comme un carre en dimension 2 ou un cube 
en dimension 3. En changeant les domaines sur lesquels ces objects sont definis, 
par exemple en considerant des vecteurs sur la sphere en identifiant chaque point 
antipodal comme equivalent, les operateurs agissant sur ces objets ont des « com- 
portements » differents. Le spectre des operateurs est entre autres different d'une 
geometrie a une autre. La geometrie spectrale etudie le spectre des operateurs 
agissant sur des objets definis sur des geometries courbes comme la sphere ou 
le tore. Comme il n'est pas necessaire de faire intervenir l'ensemble des verites 
admises concernant un enonce, la phrase precedente est certainement vraie mais 
non complete. 

Plusieurs problemes peuvent etre etudies en theorie spectrale. Ces problemes va- 
rient selon que les operateurs possedent certaines proprites comme entre autres 
etre compact, lineaire, auto-adjoint, auto-adjoint non-borne, etc. Pour certains 
types d'operateurs, comme le Laplacien A qui sera d'interet dans ce qui suit, il 
est possible de decomposer leur spectre au sens expose au chapitre I. De par la 
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forme du domaine, il peut etre possible de connaitre explicitement la forme des 
fonctions propres. 

De fagon generale, lorsque les mathematiciens ont une fonction sous leurs yeux, 
ils aiment connaitre ou elle s'annule. C'est la meme la chose avec les fonctions 
propres. L 'ensemble des zeros d'une fonction propre donnee est un ensemble ferme 
de la variete sur laquelle lesdites fonctions propres sont definies. Cet ensemble 
porte un nom particulier : il s'appelle le domaine nodal de la fonction propre 
donnee. Le complement de l'ensemble des zeros est ouvert et « decoupe » la variete 
en composantes connexes. 

Egalement, lorsque l'ensemble des fonctions propres forment une base separable 
et complete pour l'espace L 2 sur le domaine en question, il est interessant de 
s'interesser aux problemes asymptotiques. Un de ces problemes est relie a la dis- 
tribution empirique des valeurs propres. Les valeurs propres des operateurs dermis 
positifs ont des valeurs propres positives. Par consequent, etant donne une valeur 
fixe positive, combien de valeurs propres sont plus petites que cette valeur fixee 
a priori en tenant compte de la multiplicite ? 

Dans ce present manuscrit, l'analyse spectrale se fera sur le disque D de rayon 
1 ou sur un secteur S(a) d'angle a G (0, 2n) de rayon 1. dD ou dS(a) denotent 
la frontiere de D ou S(a). L'operateur d'interet sera evidemment le Laplacien 
denote A ou V 2 . Connaitre le spectre du Laplacien permet de connaitre dans ce 
cas, et dans bien d'autres cas, le spectre de l'operateur des ondes, de la chaleur 
et Schrodinger. L'espace sur lequel A est defini se denote et il s'appelle l'es- 
pace de Sobolev. Les fonctions propres formant une base pour Hq(D) avec des 
conditions de Dirichlet aux frontieres sont les fonctions de Bessel de premier type 
d'ordre 0, 1, 2, . . .. Les fonctions propres formant une base pour Hq (S(a)) sont les 
fonctions de Bessel de premier type 1,2, 

Deux problemes concernant l'operateur A sur D seront etudies plus particuliere- 
ment. Le l er probleme consistera a etudier la configuration des ensembles nodaux 
des m premieres fonctions propres Uj pour j — 1, . . . , m. Ce probleme sera equi- 
valent a calculer les m premieres valeurs propres Xj pour j — 1, . . . , m. Bref, c'est 
en voulant obtenir un algorithme efncace pour obtenir les m premieres valeurs 
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propres qu'il sera evident comment obtenir la configuration des ensembles nodaux. 
Un theoreme du a Courant donnant une borne sur le nombre de composantes 
connexes de la m e valeur propre sera l'outil principal afin d'etablir les configu- 
rations. Les ensembles nodaux seront des lignes angulaires partant du centre du 
cercle et des series de cercles concentriques. De meme, un autre probleme relie 
a l'identification de la premiere ligne nodale de la deuxieme fonction propre sur 
un secteur S(a) sera etudie. On mentionnera la valeur critique de l'angle a dans 
le cas du secteur S(a) ou la ligne decrite par l'equation en coordonnee polaire 
r = const, est preferee a la ligne decrite par l'equation 9 = const.. Nous verrons 
que si l'angle du secteur est critique, alors la l re ligne nodale n'est pas definie. De 
meme nous verrons que la multiplicite de la 2 e valeur propre A2 est double pour 
ce secteur d'angle critique. 

Un autre probleme etudie est relie a la distribution asymptotique des valeurs 
propres de A. En d'autres termes, soit Sp D (— A) = {Xj}j>o l e spectre de l'ope- 
rateur A sur le disque D de rayon 1 avec les valeurs propres ordonnees naturelle- 
ment c'est-a-dire < Ai < A2 < • • • Soit N(X), la fonction de compte des valeurs 
propres A-,- plus petites que A c'est-a-dire soit 

N(X) = £ 1 

A,<A 

Pour des valeurs de A tres gran des, 

N(X) = j-^ + R(X) (Loi de Weyl) 

ou -R(A) est le terme d'erreur. II sera d'abord reexplique clairement comment il 
a ete demontre dans [KF] que -R(A) = O(X^). Enfin, un deuxieme algorithme 
permettant d'evaluer iV(A) exactement et efficacement est donne. Cet algorithme 
n'utilise que la monotonicite des zeros des fonctions de Bessel. Comme il y sera 
explique, l'algorithme fait une « marche » qui consiste a faire des retours en arriere 
et des montees en alternance sur des paires d'entiers bien determinees par le 
probleme en soit. Les coordonnees de ces paires sont l'ordre et l'index des zeros 
des fonctions de Bessel de premier type pouvant etre les candidats possibles pour 
les valeurs propres. 



Chapitre 1 



PRELIMIN AIRES 



Voici le plan pour ce chapitre. Les references sont donnees au debut de chaque 
section. 

Dans la l re section, certaines proprietes d'interet pour ce memoire de l'operateur 
Laplacien denote A sur un ouvert regulier borne Q de M. n sont mentionnees sans 
preuve. Plusieurs definitions sont donnees avant d'y arriver permettant de bien 
definir le domaine de A et d'arriver a quelques theoremes interessants. La loi 
de Weyl (fonction de compte des valeurs propres) sera egalement presentee avec 
quelques theoremes pour l'approximer asymptotiquement sur des domaines Q 
arbitraires. 

Dans la 2 e section, il sera montre que resoudre grace a la separation de variables 
le probleme a valeurs propres Au + Xu = sur un domaine circulaire c'est-a-dire 
lorsque Q est un secteur ou un disque « engendre » l'equation de Bessel. II sera 
montre quelles valeurs A peut prendre engendrant ainsi la suite spectrale ou le 
spectre. II sera montre que les valeurs de A sont des carres des zeros de fonctions 
de Bessel. 

Dans la 3 e section, les zeros des fonctions de Bessel seront appronfondis. Plus spe- 
cifiquement, le developpement asymptotique du a Olver des zeros sera presente 
sans preuve. Ce developpement est important pour quiconque veut calculer les va- 
leurs propres de A sans avoir a resoudre les equations differentielles directement. 
Au chapitre 2, il sera montre comment en effet il est possible d'obtenir le spectre, 
les m premieres valeurs propres pour m tres grand, du Laplacien. Cette section 
introduira succinctement les fonctions d'Airy comme des solutions d'une equation 
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differentielle du meme nom. Les zeros de ces dernieres et leurs relations aux fonc- 
tions de Bessel seront abordes brievement pour nous permettre de comprendre le 
developpement asymptotique d'Olver. 

Dans la 4 e section, les ensembles des zeros des fonctions propres appeles commu- 
nement les ensembles nodaux des fonctions propres seront etudies. Un theoreme 
du a Courant sur le nombre de domaines nodaux sera enonce sans preuve. Cette 
section avec la precede nte seront utiles pour le chapitre 2. 

1.1. RAPPELS - DEFINITION DU LAPLACIEN, DECOMPOSITION DU 
SPECTRE ET FONCTION DE COMPTE 

Le lecteur retrouvera les definitions suivantes et certains theoremes sur les fonc- 
tions et les valeurs propres du Laplacien dans l'ouvrage de references [RaRa] et 
[Wi]. 

Soit Q, un ouvert borne regulier de R n . Soit C°°(fi) l'espace des fonctions definies 
sur Q differentiables indefmiment. 

Definition 1.1.1 (Le Laplacien A). Soit u G C°°(fi). Alors le Laplacien est 
I'operateur A : C°°(Q) -> C°°(fi) tel que 

. d 2 u d 2 u 

Definition 1.1.2 (Fonction propre et valeur propre). Une fonction propre 

relative a la valeur propre A est une solution non nulle du probleme 

—Am = Aw dans Q 
m = sur dfl 

Le theoreme suivant (page 115 de [Wi]) montre que les fonctions propres du 
Laplacien constituent une base de L 2 (Q). 

Theoreme 1.1.1 (Decomposition du spectre). // existe une suite non bornee de 
valeurs propres de 

< Ai < A 2 < . . . 

et une suite de fonctions propres relatives Qui est une base hilbertienne de 

L 2 (fi). 
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Les valeurs propres peuvent etre exprimees comme dans la proposition suivante 
due Rayleigh (page 115 de [Wi]). 

Proposition 1.1.1 (Quotient de Rayleigh). Pour tout n > 1, 

X n = min < / | Vu\ 2 dx | / u 2 dx = 1, / ucfjdx = pour j — 1, . . . , n — 1 > 
Jn Jn Jn ' 

Pour les resultats qui suivent concernant la fonction de compte des valeurs propres 
iVn(A), ils se retrouvent dans [Be]. 

Soit maintenant Nq(X) la fonction de compte des valeurs propres c'est-a-dire 

Nn(X) = Yl 1- 

A,<A 

La fonction de compte est une fonction a valeurs entieres definie sur les reels 
c'est-a-dire 

N n : R -> N 

qui est egalement continue a droite c'est-a-dire pour h > 

N n {\) = lim N n (\ + h) 
Si A G R est une valeur propre, alors sa multiplicite est donnee par 

lim (N n {X + h)-N Q {X-h)). 
Si A G R n'est pas une valeur propre, alors 

lim (N n (X + h)-N n (X-h)) =0. 

h->0+ 

Par consequent, si A est une valeur propre, alors il existe jsN tel que A = Xj et 
tel que 

lim N n (X -h)= lim N a (Xj - h) < j < lim N n (X + h) = lim N n (Xj + h). 

La connaissance de iVn(A) permet done de retrouver les valeurs propres Xj et 
vice- versa. 

Plusieurs resultats pour approximer N(X) ont ete obtenus. Voici certains d'entre 
eux. Le premier (voir p. 89 de [Be]) est du Weyl en 1911 qui a montre que 
Theoreme 1.1.2 (Weyl). 

Nq(X) ~ n fi(Vl)X^ lorsque X — > oo 

{27T) n 
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ou V n est le volume de la sphere de rayon 1 dans IR n et //(fi) la mesure de Lebesgue 
du domaine borne regulier Q. 

En 1952, 1953 et 1955 Levitan ecrivait une serie d'articles ou il demontrait le theo- 
reme suivant donnant une estimation du terme d'erreur. Les articles de Levitan 
sont respectivement [Lei], [Le2] et [Le3]. En 1956, Avakumovic montrait egale- 
ment le theoreme suivant dans son article [Av]. Hormander en donna egalement 
une autre preuve en 1968 dans son article [Ho]. 
Theoreme 1.1.3 (Levitan). 

^(A) = ^M^)A?+0(A^) 

Lorsque Weyl a fait la conjecture suivante qui est vraie pour tout domaine Q 
convexe et qui a ete prouvee par Ivrii dans [Iv] et, independamment par Melrose, 
dans [Melr] pour d'autres varietes riemanniennes satisfaisant des hypotheses plus 
generales. 
Theoreme 1.1.4. 

Remarque 1.1.1 (Au + \u = ou Aw + A 2 « = 0). Si au lieu d'etudier V equation 
Au + Xu = 0, V equation Au + X 2 u = est etudiee, alors il faut substituer A par 
A 2 dans les resultats qui precedent. 

Comme le disque sera l'objet geometrique d'interet dans ce memoire, si Q = D 
ou D le disque de rayon r, alors, dans [KF], il est montre le theoreme suivant 
specialise au cas du disque. Notons que [KF] etudie le probleme Au+X 2 u = 
plutot que le probleme Am + Am = 0. Le theoreme est cite presque tel 
qu'il est dans Particle arm de faciliter son explication au chapitre 3.2 
qui lui est consacre entierement. Nous comprenons bien ici que k = X 
dans le theoreme suivant. 

Theoreme 1.1.5 (Kuznetsov et Fedosov, cas du disque D de rayon r). Soit le 
probleme a valeur propre —Au = k 2 u avec u\qe> = 0. Soit la fonction de compte 
N D (k) des valeurs {k n }™ =1 ordonnees en ordre croissant c'est-a-dire soit 

N D (k) = £ 1 
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alors 

N D ( k ) = ^k 2 -^k + 0(kl) 
ou S = ixr 2 et ou L = 27ir. 

1.2. Separation de variables, equation differentielle de 
Bessel, fonctions et valeurs propres du Laplacien 
sur un domaine circulaire 

Soit S(a), le secteur du disque de rayon 1 d'angle a G (0, 2tt) ou soit D le disque 
de rayon 1. Cette section montre comment la separation de variables lors de la 
resolution de l'equation Au + \u = sur S(a) ou sur D « engendre » l'equation 
differentielle de Bessel. Les fonctions de Bessel qui sont les solutions de l'equation 
en question representent la partie radiale des fonctions propres du Laplacien sur 
un domaine circulaire. II sera montre que dans le cas du disque D les fonctions de 
Bessel impliquees dans la resolution du probleme a valeurs propres Am + Am = 
ont des ordres entiers. De meme pour le secteur S(a) il sera montre que les 
ordres des fonctions de Bessel impliques ont des ordres reels de la forme — pour 
n = 1,2, ... 

Pour un expose plus detaille de ce qui va suivre dans cette partie, le lecteur peut 
se referer entres autres a [HiCo]. En ce qui concerne les fonctions de Bessel plus 
specifiquement, le lecteur peut lire [Wa]. 

Le cas du secteur S(a) de rayon 1 et d'angle a G (0, 2ir) est etudie plus en 
profondeur dans les lignes qui suivent. Le cas du disque est presque similaire a 
celui du secteur. Soit done un secteur S(a) du disque de rayon 1 et d'angle a. Soit 
l'equation (A + X)u = avec la condition de Dirichlet aux frontieres e'est-a-dire 
u(dS(a)) = 0. En coordonnees polaires, l'equation est 

, , , „ s d 2 u Idu 1 d 2 u 

(A + xHr,e) = — + -- + - w + Xu = o 

En utilisant la separation de variables avec u(r, 9) = R(r)M(9) et en substituant, 

„2 / d 2 R , 1 dR I \ pi 

— ^ = -fr = constante = c 

R M 
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II faut done resoudre 

M"{9) + cM(6) = 
ou la solution de l'equation precedente est 

M{9) = A c cos(^c9) + B c sin(^/c9) 

Comme il faut que la fonction M soit periodique e'est-a-dire M(0) = M(a), alors 

Vc = — :=: y/c n 
a 

Puisque c ne depend que de n, il est commode de denoter par M Cn la solution de 
l'equation correspondante. Generalement, 

M Cn (9) = A Cn cosi^O) + B Cn sin(y^) , 

mais 

M(0) = M(a) 
A Cn l + B Cn = A Cn cos(nn) + B Cn 

A c = 

et, par consequent, 

M Cn (9) = B Cn sm(^9) 

Maintenant, en fixant la valeur n et en denotant par R(r) = R n (r) la solution de 
l'equation 



\ dr A r dr " "I 



tin 

II est demontre que l'equation precedente possede des solutions que pour un 
nombre denombrable de valeurs de A n (voir [HiCo]). En denotant ces valeurs 
par {A ni / C }^ 1 et leurs solutions correspondantes par R n ^(r), alors 



I dr 2 ~T r dr A n,kn n ,k J 



R 



n.k 
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Remarque 1.2.1 (Fonction de Bessel et zeros des fonctions de Bessel). La solu- 
tion k porte le nom de fonction de Bessel de premier type d'ordre — . L'ordre 
sera evident par ce qui suit. Les valeurs {X n ,k}^=i seront les valeurs propres de 
V equation originale Au + Xu = avec les conditions de Dirichlet aux frontieres. 
Les valeurs {X n ,k}kLi correspondent aux carves des zeros des fonctions de Bes- 
sel qui sont en quantite denombrable. Les fonctions de Bessel sont analytiques, ce 
qui a du sens puisque toutes les fonctions propres du Laplacien doivent etre au 
moins lisses. 

Pour plus d'information sur les fonctions de Bessel, nous pouvons lire [Wa]. 
De fagon generale, la fonction de Bessel de premier type est analytique. C'est la 
solution en serie de l'equation differentielle portant son nom et elle possede le 
developpement que voici 



Lorsque n est remplace v E C, alors J„(z) est definie par la serie suivante 



Soit done J u (z) la fonction de Bessel de premier type d'ordre reel v (l'ordre comme 
l'argument peuvent etre complexes de fagon generale) et d'argument reel z. Soit 
egalement le k e zeros de J„(z) note par jkiy). Ainsi, 





pour n G Z et zGC 




pour i/GC et zGC 






11 



Pour verifier que u n ^(r, 6) est bien la bonne solution, il est commode pour l'instant 
seulement de reecrire 



nit 
a 



A 



:=: v . 



Voici deux faits concernant les fonctions de Bessel de premier type : 



Jv+2{z) 



V 



J v (z) - J u+ i(z) 



2(z/+ 1) 



(1.2.1) 
(1.2.2) 



Par consequent, 



(J„(VAr))' = f v {V\r)(VXr)' 
= y/XJ'^y/Xr) 

= y/\(-^-J v (y/\r) - J v+1 (y/Xr) 



V 



J u (V\r) - y/\J v+1 (V\r). 



Aussi, 



(MVXr))" = 



J v {V\r) - V\J u+1 (V\r) 
rJ' v {V\r) - J u {V\r 



v 



\(j u+1 (VXr))' 



\r r 
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ce qui donne 

(A + A)J,(v / Ar)sin(^) = sm(u9) - ^ + J„(V\r) 

-V\sm(is9) (- + — + -) J„ +1 {V\r) 

\r r r J 

v 2 

-—sm(is9)J u (V\r) 
+Xsm(u9)J l/ (V\r) 



qui implique que 



(A + X)MVXr)mn(u6) = - £ + £ - £ + a) J.(VAr) sin^) 

-v^f- + — + J v+1 (y/Xr) sm(u9) 
\r r r J 

+\J l/+2 {V\r)sm(v6) 

(V 2 V V V 2 \ , i— s 

-VaY- + — + -) J v+1 (y/Xr) sm(u9) 
\r r r J 

+A^ 2( ^l 1) J i/+1 (v / Ar) - J u (V\r)^ sin(i/0) 



qui implique que 



(2 2 \ 

^-^ + ^- ^ + A- AjJ^VAr) sin(^) 



/2i/VA 2^ \f\v 
\ r r r 



^+i(v / Ar) sm(u9) 



et ainsi 



(A + A)J i ,(v / Ar)sin(z/^) = 
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Remarque 1.2.2 (Multiplicity des valeurs propres pour un secteur). Dans ce cas, 
toutes les valeurs propres pour n > 1 et k > 1 du Laplacien sur un secteur du 
disque sont simples en general sauf pour des angles critiques (voir entre autres la 
remarque 2.1.2). Les fonctions de Bessel d'ordre ne sont jamais des fonctions 
propres du Laplacien sur un secteur. 

Remarque 1.2.3 (Multiplicity des valeurs propres pour le disque). Soit D, le 
disque de rayon 1 dans IR 2 . Les fonctions propres sont encore des fonctions de 
Bessel. Cependant les fonctions de Bessel d'ordre sont admises. De fagon gene- 
rale, 

u n ,k(r,6) = J n {\/K,kr) sm(n6) ou 

u n ,k(r,0) = Jn(\/An,fc^) cos(n6) pour n > 1 et k > 1 

La multiplicity du cas correspondant a I'ordre n = est ainsi simple et celles 
correspondant aux ordres n > 1 sont doubles. 

1.3. Developpement asymptotique d'Olver des zeros des 
fonctions de bessel 

Les zeros de fonctions de Bessel sont les valeurs propres du Laplacien A sur un 
domaine circulaire, le disque ou un secteur du disque, comme il a ete discute 
precedemment. 

Plusieurs developpements asymtotiques des zeros des fonctions de Bessel et des 
zeros des derivees des fonctions de Bessel existent. Par exemples, les developpe- 
ments d'Olver, de McMahon ou de Meissel (pour J n (n)) pour ne citer que ces 
developpements les plus connus. Ces developpements sont utiles selon que I'ordre 
ou l'argument sont dans une certaine region du plan complexe. lis existent par 
exemple des developpement de Taylor pour les zeros de fonctions de Bessel d'ordre 
pour k = 0, 1, 2, . . . ou simplement pour les ordres entiers. Le lecteur peut 
lire [Ro] en ce qui concerne les developpements de Taylor relies aux zeros des 
fonctions de Bessel ou une liste de problemes ouverts est egalement presentee. 
Le developpement qui est sans doute le plus utile est celui fourni par F.W.J. 
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Olver durant les annees 50. Le developpement asymptotique d'Olver a permis en 
autre une implementation numerique stable dans plusieurs logiciels. Tout lecteur 
interesse peut lire les nombreux articles d'Olver cites dans les references comme 
[Oil], [012], [Ol5], [018] et [FaLoOl]. 

Avant de donner le developpement asymptotique d'Olver, une breve introduction 
aux fonctions d'Airy s'impose. Le lecteur peut se referer a [018], ce dernier ou- 
vrage traite des fonctions « speciales » comme les fonctions d'Airy et de Bessel. 
Une fonction d'Airy est une fonction qui est une solution de l'equation differen- 
tielle 

d 2 w 



dt 2 



tw = 0. 



L'equation precedente a deux solutions denotees par Ai et Bi qui sont nominees 
respectivement fonction d'Airy de premier type et fonction d'Airy de deuxieme 
type. 

Les fonctions d'Airy Ai et Bi sont reliees aux fonctions de Bessel de premier type 
J u avec l'ordre v — ^ comme suit 

Ai(t) = *>0 

Bid) = ^/K-Z-lCl**) +^(1* 3 )) <>° 

Les relations precedentes sont connues sous le nom de « connection formulas ». 
Les zeros de Ai seront denotes par avec k — 1, 2, . . . et > a\ > a?, > . . . 
Seule la fonction Ai sera utile pour developper les zeros de la fonction de Bessel 
J u d'ordre v > 0. II est a noter que les zeros de Ai et de Bi sont tous negatifs. 
a/- possede egalement un developpement asymptotique d'Olver (voir [018]). II est 
possible de demonter en utilisant le developpement asymptotique que 

lim 2 = 1- 

(|7r(4fc - l)) 5 

Pour la suite, nous nous concentrerons sur les zeros de la fonction de Bessel de 
premier type J u pour v > 0. En denotant par j u (k) le k e zero de J„{z) pour z > 



15 

et en denotant par 7^ = — Ofc2~t alors 



1 



■ / \ 1 3 n _i 5 — 7? , 
i k {v) = z/ + 7 fe z/sH ^iv H —v 



479 7fc 4 + 207 Vf + 2023171- 275507^ ,, + ^ 



63000 8085000 
La preuve de ce resultat est longue et complexe. Les idees de base se retrouvent 
principalement dans les deux longs articles [Oil] et [012]. 

II est possible de constater grace entres autres au developpement ci-dessus de 
j v (k) que 

lim = 1 WkeN. (1.3.1) 

v— >oo y 

II est meme prouve pour v > dans [QuWo] que 



1 



a k I ^ ■ i x Ctk I 3 2 23 

z/-- T z/3< Jfc (z/)<z/-— 1/3+ a fc — . (1.3.2) 

23 23 /U //3 

Grace a (1.3.1) et a (1.3.2), il est possible d'evaluer l'ordre de grandeur des zeros. 
Ceci sera utile ulterieurement (voir les remarques A. 5.1 et 3.3.4) afin de connaitre 
le degre de precision requis pour evaluer numeriquement les zeros des fonctions 
de Bessel. 

1.4. Les ensembles nodaux des fonctions propres et leurs 
structures, cas du disque et du secteur en exemples 

Le theoreme principal de cette section concerne le nombre de domaines nodaux 
qu'une fonction propre peut avoir. II est enonce entre autres dans [Sh]. Avant 
d'arriver au theoreme, un rappel de ce qu'est l'ensemble nodal d'une fonction 
propre s'impose. De meme, les exemples pour le secteur du disque S(a) de rayon 
1 et d'angle a G (0, 27r) et le disque D de rayon 1 sont donnes. 
L'ensemble nodal d'une fonction propre est son ensemble de zero. Si la dimension 
du domaine est d, alors la dimension des ensemles nodaux est d — 1. En d'autres 
termes, pour la j e fonction propre (fj defini sur un domaine ouvert assez regulier 
Q C lR d , si Zj est l'ensemble des zeros, alors 



Zj = ^y?j(x) = 0, x E fij. 
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Zj est un ensemble ferme. 



Dans le cas d'un secteur S(a) de rayon 1 et d'angle a G (0, 27r), une fonction 
propre typique est donnee par 

[nir6 s 



Par consequent, 

Un,k(r,9) = 4» J^{^/\ n ,kr) = ou sin ^— — J = 

Ainsi 



"rvKt 

sin ' 

a 



t 



dn pour d = 0, 1, . . . , n 

da 
n 

lequel cas donne une suite de lignes angulaires decrites par les multiples entiers 
d = 0, 1, . . . , n de l'angle -. Le cas ou 

JiHL(A/A^r) = 
t 



pour j = 1, . . . k 



donne une suite de cercles concentriques avec le cas correspondant a la frontiere 
lorsque j = k, c'est-a-dire lorsque r = 1. 

L'ensemble des zeros d'une fonction propre decoupe le domaine de definition des 
fonctions propres en composantes connexes ou la fonction propre alterne de signe. 
Dans le cas du disque ou d'un secteur circulaire, les composantes connexes ont la 
forme de secteurs circulaires comme il est possible de le voir sur les figures de la 
section 2.2.1. Une image vaut milles mots. 
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Remarque 1.4.1 (Nombre de composantes connexes). 
(cas d'un secteur) 

Pour une fonction propre donnee u n .k{r,6), le nombre de composantes connexes 
sur lesquelles la fonction propre alterne de signe est donnee par nk. 
(cas du disque) 

Pour les fonctions propres sur le disque de multiplicite 1 c'est-a-dire pour celles de 
la forme Uo,fc(^ ? ^)? Is nombre de composantes connexes est k. Pour les fonctions 
propres de multiplicite 2 c'est-a-dire celles de la forme u n ±(r, 9) avec n = 1, 2, . . 
le nombre de composantes connexes est 2nk. 

Avant de conclure cette section, voici un theoreme celebre du a Courant qui sera 
tres utile pour la prochaine section qui peut etre retrouve entre autres dans [Sh]. 
Theoreme 1.4.1 (Theoreme de Courant sur les ensembles nodaux). Soit Q C M. n , 
un domaine tel que dQ est suffisamment regulier par morceaux. Soit X m , la m e 
valeur propre. Alors pour une fonction propre u m correspondant a \ m , le nombre 
de domaines nodaux est au plus m. 



Chapitre 2 



LES ENSEMBLES NODAUX DES FONCTIONS 

PROPRES 

Soit D ou S(a), le disque de rayon 1 ou un secteur de rayon 1 et d'angle a G 
(0, 2tt). Soit A l'operateur differentiel de Laplace. Soit A, une valeur propre A 
telle que (A + X)u = avec u G C°° etant la fonction propre associee a A. 
Dans un premier temps, une conjecture concernant la l re ligne nodale du probleme 
de Dirichlet sur un domaine bidimensionnel est etudiee. La conjecture en question 
a ete soulevee par Payne et affirme que pour tout domaine bidimensionnel, la 
l re ligne nodale touche toujours la frontiere. Deux cas extremes s'offrent a la 
conjecture, le rectangle et le secteur circulaire. Le cas du rectangle etant tres 
simple puisqu'il est facile d'identifier ou se trouve exactement la premiere ligne 
nodale comme il sera montre plus loin. Le cas du secteur etant un peu plus 
interessant, cela motivant en soi la prochaine partie. A la fin de la f re section, 
le lecteur trouvera une proposition donnant Tangle critique a pour lequel tout 
secteur d'angle a < a a une ligne nodale en coordonnees polaires de la forme 
r = const et pour lequel tout secteur d'angle a > «o a une ligne nodale en 
coordonnees polaires de la forme 6 = const. De meme, lorsque a = ao, nous 
verrons que la ligne nodale est indeterminee. 

Dans un deuxieme temps, soit un entier m > 1 generalement tres grand comme 
m—1- 10 6 par exemple. Comment construire la suite des valeurs propres {Xj} 1 JL 1 
en ordre croissant de l'operateur A sur D ou S(a) ? En d'autres termes, il faut 
etablir un algorithme efficace pour determiner le spectre jusqu'au rang m. Par 
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efficace, cela veut dire ne pas resoudre l'equation differentielle avec des methodes 
comme les elements finis qui ne permettraient pas de construire la suite spectrale 
pour des valeurs tres grandes comme to — 1 • 10 6 et ni meme d'identifier les 
ensembles nodaux de fagon exacte. Par efRcace, il faut se servir de toute la theorie 
decrite jusqu'a maintenant et du fait que dans certains logiciels, en l'occurence 
Maple, il existe des routines tres stables permettant de trouver les zeros des 
fonctions de Bessel. La construction de la suite spectrale jusqu'a to est equivalente 
a denombrer le nombre de composantes connexes disjointes obtenues pour chaque 
Uj en enlevant du domaine original les ensembles nodaux des fonctions propres 
Uj pour j = 1,2, ... , to. Les valeurs propres sont donnees en fonction de deux 
parametres n et k entiers non-negatifs et connaitre le spectre jusqu'a l'entier to, 
consistera a etablir une bijection entre N 2 et N. 



2.1. La premiere ligne nodale de la deuxieme fonction 

PROPRE 

Dans cette partie, la premiere ligne nodale de la deuxieme fonction propre de A 
sur un domaine convexe du plan est etudiee. Cette ligne est l'ensemble des zeros de 
la deuxieme fonction propre w 2 du Laplacien d'un domaine convexe du plan, avec 
des conditions de Dirichlet sur la frontiere. Pour la premiere fonction propre, u±, 
elle ne possede qu'une seule ligne nodale qui s'annule sur la frontiere du domaine 
en question et qui, par consequent, est triviale. Quant a w 2 , par le theoreme de 
Courant sur les ensembles nodaux, elle possede au plus 2 lignes ; la premiere etant 
toujours triviale et s'annulant sur la frontiere et quant a la deuxieme ligne, non 
triviale celle-ci, elle peut toucher a la frontiere ou etre entierement contenue a 
l'interieur du domaine. 

Mathematiquement parlant, soit Q un domaine convexe du plan. Soit -u 2 , la 
deuxieme fonction propre de A, c'est-a-dire Au 2 + A 2 w 2 — et u 2 (dQ) = 0. 
Soit r = {z G Q,u 2 (z) = 0}, alors deux possibilites s'imposent : T H dQ = ou 
P n dQ j£ 0. Dans le cas ou f n dQ = 0, alors T CQ. 
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L.E. Payne a conjecture que T D dQ ^ pour tout domaine du plan. S.-T. Yau a 
conjecture que T PI dQ 7^ pour tout domaine convexe du plan. Pour un domaine 
borne convexe Q C M 2 , Melas en 1992 dans [Mela] a montre le theoreme suivant. 
Theoreme 2.1.1 (Melas). Soit Q C IR 2 un domaine borne convexe avec des fron- 
tieres CP° , alors la ligne nodale de la deuxieme fonction propre u 2 doit intersecter 
la frontiere dQ exactement en deux points. 

Pour un secteur S(ct) de rayon 1 qui est un domaine borne convexe du plan, la 
ligne nodale de la deuxieme fonction propre peut done toucher dS(a) de deux 
fagons. La ligne nodale intersecte la partie de dS(a) decrite par r = 1 ou les deux 
parties decrites par 9 = et 9 = a. 

Question. Quelle doit etre la valeur de a faisant en sorte que nous avons une 
situation plutot qu'une autre ? 

Remarque 2.1.1 (Deux possibilites pour -u 2 )- a done deux sous domaines 
simplement connexes disjoints separes par une ligne nodale, radiale ou angulaire. 
Les seuls zeros admis des fonctions de Bessel sont Ji(-) , ii(2-) etj2{-). Puisque 
^i = 3i\~) necessairement, alors les deux seules autres possibiltes sont ji(2^) 
et jaf—) chacun donnant respectivement une ligne nodale angulaire et une ligne 
nodale radiale et, dans chaque cas, deux sous domaines simplement connexes oil 
U2 alterne de signe. 

La fonction jk{v) avec v = - etant lisse et croissante avec l'ordre u, alors la ligne 
nodale radiale se produit si et seulement si 



et l'equation de la ligne radiale est donnee par 



r = 




7T 

avec v — — . 

a 



Quant a la ligne nodale angulaire, elle se produit si et seulement si 



et l'equation de la ligne angulaire est donnee par 



9 



a 



2 
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Lorsque v = - satisfait l'equation ji(2u) = j2{v), alors l'angle a est critique et 
la deuxieme valeur propre a une multiplicite double. Si nous resolvons l'equation 

j 2 (z/)-j 1 (2z/) = 0, 

alors la solution est 

v « 2.823823 a « 1.112531 « 63.743334 degrees 
et la l re ligne nodale n'est done pas definie. En effet, soit 

/ 71 \ \ . /7t0\ 



%,2 



= J*rfjif2— ) r )sinf— ) (2.1.2) 
les deux fonctions propres associees respectivement aux zeros et ii(2^). 

\ an / V ( 



Comme 

-) 

alors n'importe quelles combinaisons lineaires des fonctions propres (2.1.1) et 
(2.1.2) est une solution de l'equation A(u) + \u = sur S(an). 

Remarque 2.1.2 (l re ligne nodale et secteur critique). Etant donne un angle a et 
le secteur du disque correspondant d 'angle a, les valeurs propres sont generalement 
de multiplicite 1. Soit l'angle critique an et I'ordre associe = ~ tel que ^2(^0) — 
ji(2z/o) = 0. Lorsque l'angle a = an, la deuxieme valeur propre a une multiplicite 
de 2. La premiere ligne nodale n'est done pas definie pour cette valeur critique a 
puisque n'importe quelles combinaisons lineaires des fonctions propres (2.1.1) et 
(2.1.2) est une solution du probleme A (it) + Xu = sur S(ao). 
II est interessant de visualiser les ensembles nodaux de la 2 e valeur propre pour 
quelques valeurs a ou, en d'autres termes, la l re ligne nodale. Nous voyons en 
effet que lorsque a — >■ ckq, la ligne nodale est indeterminee, en effet, Matlab prend 
au « hasard » une combinaison lineaire des deux fonctions propres. Les angles 
mentionnes en-dessous de chaque graphique sont exactes et ne sont 
pas approximatives. Le but etant de donner a Matlab des angles se 
rapprochant de plus en plus pres de an. 
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3hfcnrdferghl.il ^91.0072 




0.1 M 03 0.4 05 OS 0.7 Oi OS 1 



Figure 2.1. l re ligne nodale, a = 1.11, A 2 « 91.0072 



Sectaurdfengle 1.112531 ^=90.8402 




0.1 Oi 03 0.4 05 OS 0.7 0* OS 1 



Figure 2.2. l re ligne nodale, a = 1.112531, A 2 « 90.8402 
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Ssotaurd'argh1.11253100a739J(t1 ^=S0.S371 

os r 





0.1 02 03 0.4 05 OS 0.7 tH OS 1 



Figure 2.3. l re ligne nodale, a = 1.112531602739001, A 2 « 90.8371 



Smteardktgfe 1J! ^61.9144 

1 r 




0.1 (B 03 0.4 05 OS 0.7 tB OS 1 

Figure 2.4. l re ligne nodale, a = 1.2, A 2 « 81.9144 



Voici done a la lumiere de ce qui precede la petite proposition suivante concernant 
l'« emplacement » de la premiere ligne nodale de la deuxieme fonction propre. 
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Proposition 2.1.1 (« Emplacement » de la premiere ligne nodale). Soit S(a) un 
secteur du disque de rayon 1 et d 'angle a. Soit cvq la solution de I 'equation 

M-H'(-) 

Vao/ \ao/ 

oiijk{v) est le k e zero de la fonction de Bessel de premier type d'ordre v > 0. Alors 
si V angle a du secteur S(a) est superieur a ao, alors la premiere ligne nodale de 

■ 2 ( \ 

la deuxieme fonction propre est de la forme 9 = | sinon elle de la forme r = 
avec v — -. 

a 

Demonstration 2.1.1. Lire la remarque 2.1.1. 

2.2. ALGORITHME POUR ORDONNER LES VALEURS PROPRES ET 
DETERMINER LA STRUCTURE DES ENSEMBLES NODAUX 

Question Supposons donne un entier m > 1. Est-il possible de determiner effi- 
cacement le spectre du Laplacien sur le disque D ou un secteur S(a) de rayon 
1 jusqu'a l'entier m? En d'autres termes, comment obtenir toutes les valeurs 
propres \j pour j — 1, . . . , m telles que 



< Ai < A 2 < . . . < A 



in 



Reponse Repondre a la question precedente, c'est connaitre la bijection entre 
j\{y) et A m ou v = n pour le disque D et oil v = — pour le secteur S(a). 
Nous avons done un probleme de classification des valeurs j1{k). Etant donne m, 
quelles valeurs de n et k sont admissibles en fonctions de m? Dans ce qui suit, 
nous ecrirons v = u(n). 

Idee de l'algorithme Par le theoreme de Courant, il faut trouver toutes les 
paires (n, k) G N 2 telles que nk < m, ce sont les paires admissibles. De meme, la 
connaissance de toutes ces paires (n, k) implique la connaissance des paires (n', k') 
telles que n'k' < m' pour m' < m. Par consequent, evaluer j„( n \(k) pour toutes 
les paires (n,k) telles que nk < m assure d'obtenir {A J }™ =1 ^ {jl( n ){k)}nk<m- H 
ne reste done qu'a ordonner et a prendre en compte la multiplicite dans le cas du 
disque. Dans le cas de S(a), il ne faut qu'ordonner lineairement {jl^{k)} n k<m 
et prendre les m premieres valeurs, ce qui donne < Ai < . . . < A m . Dans le 
cas du disque, il faut regarder les valeurs ou n = et celles n / et extraire 
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les valeurs de {jl^(k)} n k<m en repetant deux fois celles ou n ^ pour obtenir 
< Ax < ... < A m . 

L'algorithme precedent a ete implements en Maple et est donne en annexe. 
Remarque 2.2.1 (Ameliorations possibles). Comme a n'importe quelle pro- 
gramme informatique, plusieurs ameliorations peuveut etre apportees. 

1) La premiere, pour ceux dont les ressouces le permettent, seraient de paralleliser 
le code avec MP I. 

2) Afin de reduire la taille en memoire utilisee et egalement de permettre de 
sauvegarder en mode binaire (car Maple ne fait pas la distinction entre le mode 
binaire et ascii sur les stations Unix), reprogrammer en C+ + avec par exemple la 
librairie NTL (Number Theory Library) compilee avec GMP (Gnu Multiple Preci- 
sion arithmetic) procureait sans aucun doute des economies de temps et d'espace. 
Cela implique cependant d'implementer une fonction pour le calcul des zeros de 
fonctions Bessel de premier type qui serait stable, ce qui n'est pas necessairement 
evident a faire efficacement. 

Remarque 2.2.2 (Avantages de l'algorithme). A la connaissance de Vauteur, 
l'algorithme donne ici est un des plus efficace pour determiner les valeurs propres 
et les fonctions propres de I'operateur A sur le disque. II ne requiert aucunement 
de resoudre V equation differentelle Au + \u = sur le disque. Calculer les ze- 
ros des fonctions de Bessel etant beaucoup plus rapide et stable que de resoudre 
V equation differentielle, l'algorithme est d'autant plus rapide et efficace. Egale- 
ment plusieurs des methodes numeriques pour resoudre I' equation differentielle 
requiert I'utilisation de matrices done le nombre d' entrees est proportionnel au 
carre du nombre de valeurs propres desirees et, par consequent, il serait inutile 
d' essay er de calculer par exemple les 1-10 6 premieres valeurs propres en resolvant 
V equation differentielle avec ces methodes numeriques, ce qui par contre a ete fait 
raisonnablement avec l'algorithme. 

2.2.1. Resultats comparatifs 

Le spectre du disque pour les valeurs propres {Aj}^ 6 ! a ete calcule environ en 
114576.4710 secondes de calcul soit environ 31 heures 50 minutes sur un simple 
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PC a deux processeurs avec 2GB de RAM. La liste d'une taille de 30 MB est 
disponible aupres de l'auteur. Dans ce ficher texte, la l re colonne indique le rang 
de la fonction/valeur propre, la 2 e colonne indique la valeur propre, la 3 e indique 
l'ordre de la fonction de Bessel et la 4 e indique l'index du zero de fonction Bessel. 
Dans ce qui suit, pour chacune des valeurs prop res listees dans les tableaux, la 
representation en courbes de niveaux de la fonction propre correspondante est 
montree plus loin sur les graphiques. Grace aux courbes de niveaux, il est facile 
d'y compter le nombre de domaine nodaux et de comparer avec les valeurs de 
k (index) et n (ordre) dans les tableaux. Les courbes de niveaux des fonctions 
propres ont ete obtenues a l'aide de Matlab en solutionnant les equations diffe- 
rentielles a l'aide de la methode des elements finis. Par exemple, a des fins com- 
paratives, les 100 premieres valeurs propres et leurs fonctions propres associees 
ont ete obtenues en pres de 3 jours de calcul en Matlab. En plus de prendre un 
temps excessivement long, la visulation des courbes a niveau devient tres difficile 
lorsque le rang des fonctions propres est eleve. 

Des programmes en Matlab sont egalement donnes en annexe. Ces programmes 
permettent de resoudre les problemes de valeurs propres de meme que de visualiser 
les ensembles nodaux des fonctions propres. 
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m = rang 


\ m (Matlab) 


jl(n) (Maple) 


n = ordre 


k = index 


3 


149.4956 


149.4529 


2 


1 


6 


278.9782 


278.8316 


3 


1 


7 


310.5226 


310.3223 


1 


4 


16 


646.8263 


646.0310 


5 


1 


25 


991.8437 


989.7291 


3 


5 


27 


1088.4330 


1085.9440 


4 


4 


30 


1157.8175 


1155.2319 


5 


3 



Table 2.1. Resultats pour le secteur S(a) avec a = | de rayon 1 



m = rang 


\ m (Matlab) 


j 2 k (n) (Maple) 


n = ordre 


k = index 


3 


23.0006 


22.9968 


3 


1 


9 


58.4266 


58.4019 


7 


1 


11 


69.3865 


69.3521 


8 


1 


19 


108.3025 


108.2183 


2 


3 


26 


152.1173 


151.9596 


14 


1 


29 


165.6505 


165.4521 


5 


3 



Table 2.2. Resultats pour le secteur S(a) avec a = 2exp(l) ss 
5.4366 de rayon 1 

Les figures 2.5 et 2.6 representent les domaines nodaux des fonctions propres 
correspondantes aux valeurs propres du tableau, le nombre de composantes 

connexes est donne par nk. 



Sedeurdartgte0.7S54 ^=1496237 Sedeurdartgte0.7S54 7^=311 1235 

as, 

07- 



OB 
V3 

^^^^^ 

1 




0.1 02 03 0.4 QjS OB 0.7 Ofi 03 1 
Secteurd'angls 0.7854 9^=649 2317 

Ofi 
0.7 
OB 

« » 

1 


0.1 02 03 0.4 05 0B 0.7 Ofi 03 1 
Sedeurd'angle 0.7854 9^=998.1934 

A 


0.1 02 03 0.4 0*5 OB 0.7 Ofi 03 1 
Secteurd'arigle 0.7854 ^1095.91-86 

Ofi 
0.7 
OB 

Si il 


0.1 02 03 0.4 05 0B 0.7 Ofi 03 1 
Secteur da rgte 0.7854 ^=1165.5948 

▲ 

♦ » 

• •• 


0.1 02 03 0.4 05 OB 0.7 Ofi 03 1 0.1 02 03 0.4 05 OB 0.7 Ofi 03 1 

Figure 2.5. Domaines nodaux de la 3 e , 7 e , 16 e , 25 e , 27 e et 30 
fonction propre pour le secteur S(a) avec a = f de rayon 1 
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S«teJ[d'arBte5,4366 ^=23.0122 

1r 1 

obL OB 
0.4 1 0.4 
Oli j 02 

02 1 -02 
0.4 L -0.4 
0.6 L -0.6 

o.eL -o.e 
_il 1 1 1 1 1 1 1 1 1 _i 


Sffitetird'argts 5.4366 ^=53.5007 


-1 -0.fi -O.S -0.4 -0.2 02 0.4 OB OS 1 
Secteurd'ari3te5.43S5 ^,=69.4396 

'[ gm ^ 

. if I 

0.4 1 0.4 
02L ^^H^k 

J 

0.2 1 -0.2 

0.4 \ -0.4 

o.e I -o.e 
_il 1 1 1 1 1 1 1 1 1 _i 


-0.fi -O.S -0.4 -0.2 02 0.4 0B OS 1 
Secteurd'an3b5.43ffi ^=106^544 

* 


-1 -0.fi -O.S -0.4 -0.2 02 0.4 0B 0.8 1 

1r ' 1 

0^L 05 

0.4 L 0.4 
02L 1 02 

'A 

0.2! ^ -0.2 
0.4 L -0.4 

_-|l 1 1 1 1 1 1 1 1 1 _1 


-0.fi -O.S -0.4 -0.2 02 0.4 0B OS 1 
Secteurd'an3b5.43ffi ^=1662459 


-1 -0.fi -O.S -0.4 -0.2 02 0.4 0B OS 1 

Figure 2.6. Domaines nodaux de h 
fonction propre pour le secteur 5(a) a 


-0.fi -O.S -0.4 -0.2 02 0.4 0B OS 1 

i 5 e , 9 e , ll e , 19 e , 26 e et 29 
vec a ~ 5.4366 de rayon 1 



30 



m = rang 


jl(n) (Maple) 


n = ordre 


k = index 


1 


5.783186 





1 


4 


26.374616 


2 


1 


5 


26.374616 


2 


1 


6 


30.471262 





2 


7 


40.706466 


3 


1 


8 


40.706466 


3 


1 


28 


135.020709 


2 


3 


29 


135.020709 


2 


3 


30 


139.040284 





4 


31 


149.452881 


8 


1 


32 


149.452881 


8 


1 


33 


152.241154 


5 


2 


34 


152.241154 


5 


2 










1000000 


4004017.840283 


1533 


69 


2000000 


8005695.299643 


391 


714 


3000000 


12006894.927313 


305 


955 


4000000 


16008031.629743 


1820 


498 


5000000 


20008771.752617 


1624 


707 


6000000 


24009692.370612 


713 


1220 


7000000 


28010694.266301 


3663 


276 


7912680 


31685063.008767 





1792 



TABLE 2.3. Resultats pour le disque 

Les figures 2.7 et 2.8 representent les domaines nodaux des fonctions propres 
correspondantes aux valeurs propres du tableau, le nombre de composantes 
connexes est donne par 2nk si n > 1 (multiplicite 2) et par k si n = 

(multiplicite 1). 



Deqje ^=5.7843 Deque 




1 1 1 1 1 1 1 1 1 1 1 -l' 1 1 1 1 1 1 1 1 1 

-1 -0.£ -0.6 -0.4 -0.2 02 0.4 OB OS 1 -1 -M -0.S -0.4 -0.2 02 0.4 OB OS 1 



Figure 2.7. Domaines nodaux de la l re , 4 e , 5 e , 6 e , 7 e et 8 e fonction 
propre pour le disque D de rayon 1 



Chapitre 3 



LOI DE WEYL POUR LE DISQUE 



Dans cette section, une explication approfondie de Particle [KF] est donnee. 
C'est dans cet article qu'il a ete montre comment la loi de Weyl pour le disque D 
de rayon r est dans l'ordre 2/3. Je tiens a remercier Igor Wigman et Dominique 
Rabet pour les discussions tres enrichissantes que nous avons tenues afin d'eclaicir 
les multiples points tres obscures de Particle. 

L'article montre que la loi de Weyl, en debutant avec le probleme Am + k 2 u = 0, 
pour le disque D de rayon r est dans Pordre 2/3. En d'autres termes, soit le 
probleme a valeur propre —Am = k 2 u avec u\qd = 0. Soit la fonction de compte, 
Nu(k), des valeurs propres {k n }^ =1 ordonnees en ordre croissant c'est-a-dire soit 

N D (k) = 1 

alors il est montre que 

N D ( k ) = ^k 2 -^k + 0(kl) 

47T 47T 

ou S = nr 2 et ou L = 2itr. 

Afin d'alleger la notation, nous ecrirons simplement N(X) au lieu de 
Nd(X). Seulement pour l'explication du theoreme de Kuznetsov et Fe- 
dosov, nous utiliserons k au lieu de A afin de conserver la notation 
identique a celle de l'article facilitant ainsi les references. 

Avant de commencer a expliquer Particle, je tiens a faire un « detour » theorique 
sur Papproximation uniforme de la fonction de Bessel de premier type, J n (x) pour 
x > n > 0. Les fonctions de Bessel de premier type d'ordre n = 0, 1, 2, . . . sont au 
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coeur de l'etude du spectre du Laplacien sur le disque. Ce « detour » se justifie 
d'une part parce que les auteurs donne une formule asymptotique pour J n (x) pour 
laquelle il faut travailler quelque peu avant d'y arriver et, d'autre part, parce que 
la formule donnee n'est pas effective tandis que je montre comment obtenir une 
expansion uniforme asymptotique effective c'est-a-dire sans terme O. La theorie 
sur les expansions uniformes des fonctions a ete principalement developpee par 
F.W.J. Olver dans les annees 1950. 

A la fin de ce chapitre, je donne un algorithme permettant de calculer emcacement 
et exactement N(X). En utilisant la monotonicite des zeros des fonctions de Bessel 
ainsi qu'un principe de « marche » qui consiste a faire des retours en arriere et des 
montees en alternance sur des paires d'entiers bien determines par le probleme, 
nous verrons qu'il est possible d'evaluer N(X) emcacement sans avoir a calculer 
nai'vement toutes les valeurs propres telles que Xj < A. Par methode naive, je 
veux signifier Futilisation de l'algorithme de la section 2.2 afin de determiner la 
valeur du rang maximale de la valeur propre correspondant a la valeur de iV(A). 

3.1. Approximation uniforme de la fonction de Bessel de 
premier type 

Le lemme suivant decoule de deux theoremes qui seront exposes bientot. La par- 
tie la plus difficile a ete d'evaluer les variations totales de certaines fonctions qui 
apparaitront bientot. II est a noter que les auteurs de [LaWo] affirment avoir 
analytiquement trouve les points stationnaires des fonctions en question que vous 
pouvez consulter a la section 3.1.1 sans le montrer dans leur article, ils ne se 
contentent que de donner les valeurs numeriques de variations totales qui sont 
evidemment les memes que les miennes sans meme donner explicitement les fonc- 
tions comme je l'ai fait. Certaines des fonctions impliquees permettant de borner 
l'approximation d'Olver avaient deja ete obtenues par Olver, je me suis contente 
de les recalculer numeriquement pour obtenir plus de decimales qu' Olver a pu le 
faire avec les ordinateurs des annees 50. Les resultats utiles pour cette section se 
retrouvent dans [013], [014] et [016]. 
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Lemme 3.1.1. Soit x > n > 0, J n (x) la fonction de Bessel de premier type et 



7r - 2 2^ 



f n (x) = (l + 6 3 ) x /-(x 2 -n')*J n (x) (3.1.1) 



alors I'inegalite suivante est satisfaite 



f n {x) - cos \ 7) - 



Ui So 
<— + - 

77 n 



i | n | S 1 1HS! \R 2 \s 1 \vi\s 2 \R' 2 \S2 . , , 

+ \R 2 \ + — j H ^— H 5 1 1 h V7r e 3 

rr 77/77 77/ 7177 n 

oil s%, s 2 , u%, v 1, R 2 , R' 2 , 5 3 , e 3 soni exposes plus loin et ou rj est comme suit 

r] = Vx 2 — n 2 — n arccos ^— ^ 
De meme, pour une valeur fixe arbitraire de x, il existe n x < x telle que 



f n (x) = cos (77 - Y) 



01 1W 1 



.71/ V77, 

Remarque 3.1.1 (n x ). La valeur n x peut etre bien approxime en utilisant Vex- 



cellente approximation de fj(y) = y 1 — y 2 — y arccos(?/) pour y G (0, 1) qu'est 

(l-y)l 3 r- 

En effet, de par I'inegalite de la remarque precedente, qui m'a ete fournie par Igor 
Wigman, avec < y — - < 1, 

(1-1/) 5 < ^77(7/) ^ < ^(~) = ^("'^ 

et, par consequent, 

Ax / n\f , . /Sv^xi , n /3\/2\f 1 

K^il--) <^,n)^(— ) <l--^n<x-(— ) *a 

II est a noter que x — x% < x — (^p) l x^ pour 2 > (^p) 6 ~ 1.423828. 
De fagon equivalente comme il en ressortira du theoreme 3.1.1, 

fn(x) = e 1 (r l ) + e 1 (r])^p- + -e 2 (r ] )(-O^B (0 + ^ 3 (3.1.2) 



n 2 n 



ou 



#1 (77) = cos ^77 - ^7rj + sin ^77 - ^7rj — |- + -R 2 



4 / rj 1 



6 2 {r)) = sin (77 - -7r) - cos (77 - -7r) p- + # 2 
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alors 

Jn{x) = TT^Xl^ ~ n 2 )~^fn{x) 

1 + S V 7T 

Remarque 3.1.2 (Denombrement des zeros de J n {x)). Seule la connaissance de 
f n {x) est necessaire au denombrement des zeros de J n (x). 

Avant de debuter la preuve, voici les deux theoremes qui seront utilises pour 
prouver le lemme et se retrouvant respectivement dans [017] et [FaLoOl]. 
Avant d'enoncer les theoremes, voici les definitions des fonctions qui apparaitront 
dans les theoremes. Dans ce qui suit, Ai(x) et Bi(x) denotent respectivement la 
fonction d'Airy de premier type et deuxieme type. Ainsi, 

„, N fBi(x)Y 
E[x) = - , . pour c < x < 00 

\Ai(x)J P " " 

E(x) = 1 pour —00 < x < c 

M(x) = (E 2 (x)Ai 2 (x) + E~ 2 (x)Bi 2 (x))^ pour iGl 

A = sup(7r|x|^M 2 (a;)) 

w 1.039522542988 a x w 1.321915092767 



ou c est la plus grande valeur negative telle que Ai(x) = Bi(x). Nous nous 
entendons bien que A n'a aucun lien avec les valeurs propres. C'est la 
lettre utilisee par Olver pour designer le supremum ci-dessus. 

Quant au symbole V[ a ^){g)i il denote la variation totale d'une fonction g sur un 
intervalle (a, b) C dom.{g). Pour plus d'information, voir le l er chapitre de [018]. 
Theoreme 3.1.1 (Approximation uniforme de J n {nz)). Soit J n {nz) la fonction 
de Bessel de premier type d'ordre n e ffi. et d 'argument nz avec z G C, alors 
V expansion uniforme de J n (nz) est donne par 

■«») - tt^Mt^)\^4^ {3A3) 
+ ^Ef +£ J (3,.4) 
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oil les termes d'erreur 821+1 et €21+1 sont donnes comme suit 




Le lecteur interesse a connaitre les recurrences permettant de calculer 
A s et B s dans le theoreme qui suit peuvent consulter [014]. J'ai trouve 
Ai(Q et Bi(Q en resolvant les recurrences, mais cela ne requiert qu'une 
bonne dose de patience et je me suis contente de les donner a la section 
3.1.1. ( est donne plus loin a l'equation (3.1.16). 

Remarque 3.1.3 (n et z). Dans ce qui suit, seuls les cas ou n G N et z G M + 
seront utiles. 

Theoreme 3.1.2 (Approximations uniformes de Ai et Ai'.). Soit Ai la fonction 
de Airy de premier type sa derive Ai' . Soit x > et = £. Alors 




(3.1.7) 




§s-ij 



(3.1.8) 



ou 




et 




(3.1.9) 




(3.1.10) 



ou 
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avec 



2T(3s + i 

u s = ■= ^— 2 —- s>0 3.1.11 

(6s + 1) , 

■ - J -u s s > (3.1.12) 



(6s 

M = 1 (3.1.13) 
Voici maintenant la preuve du lemme 3.1.1 

Preuve du Lemme 3.1.1. A partir des equations (3.1.3), (3.1.5) et (3.1.6) du 
theoreme 3.1.1 avec I = 1 dans 1' approximation de J n (nz), nous avons que 

1 1 / 4C / 2 ,( M0\ Ai'(n*Q 



avec 



0.013227 /0.218423\ 
63 < 5 exp 3.1.14 

n 6 \ n J 

0.008352 /0.218423\ 
M < = exp . 3.1.15 

n 6 \ n J 

Nous reviendrons a la section 3.1.1 sur les evaluations numeriques de 63 et 63. 
Maintenant, un changement de variable qui « defait » l'homothetie nz est donne 
par x = nz. En se concentrant seulement sur z > 1 c'est-a-dire sur z G (l,oo), 
alors (voir [016]) 

2 

3 \ 3 2 1 

rfi— x>n>0 (3.1.17) 

n,3 



Oil 



n 



Ainsi 



V x 2 — n 2 — narccos ^— j x > n > 

Ai(nk) = Ai(-(nf(-C))) 
Ai'(nH) = ^'(-(nl(-C))) 
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En substituant S = 2 dans les approximations de Ai et Ail donnees respective- 
ment par (3.1.7) et (3.1.9), et en denotant pour des raisons de commodite Q\{j]) 
et 9 2 (t]) par 

0i (77) = cos (p - ^tt) + sin (r) - -^j ^ + R 2 (3.1.18) 
6 2 {rj) = sin (t;-^)- cos (77-^)^ + ^ (3.1.19) 

* eX K4)w (3 - L20) 
M * (3 - L21) 



alors 



^Mc) = -^—-—^1(77) 

Ai'(n*0 = -pnB(-C)^ 2 (»7) 

\/7T 



Par consequent, pour 2 > 1 ou bien £ < 



Unz) = T ^J r ( r ^i_y( r _LJ r e 1 (T / )(l + 



?>, 2 



+ 4=^ a (i7)(-o^o(c)+c3 

V 71 " r?- 3 



et done pour x > n > avec £ < 0, 



1 + 03 V vr V n z 



Maintenant soit (voir [016]) 



Sl = sup 1^(01 = -^- (3.1.22) 

C<o 

s 2 = |C(2i)|5|5 (2i)| (3.1.23) 

« 0.010862854400 (3.1.24) 

et ou zi est le seul point stationnaire de \((z)\i Bo(((z)) pour z > 1 £ < tel 
que 

21 = 1.979495483061 avec COi) = -1.000459796360 . 
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Par consequent, 



\Bi(v)Ai{C)\ < \0M\si 



n 2 n 2 ^ n 2 



-|^ 2 (r7)(-C)^o(C)l < -|fcfa)|a a 



< g 2 | \vi\S2 | |-R 2 l g 2 

~ n nrj n 



En denotant par 



2 



n 2 n 



alors pour x > n > 0, l'inegalite suivante est satisfaite 



f n (x) - COS (p-J 



rj n 



1 f \ t> \ 1 s i , , l-Raki , M s 2 . |^ 2 I S 2 . rx 

+ i? 2 + — j H 5— H 5 1 1 V V7T e 3 

\ rr 72/77 77/ tit? n 

Ceci termine la preuve de la premiere inegalite. Pour la deuxieme partie 

du lemme, procedons comme ce qui suit. 

Soit la fonction f n (x) donnee par l'equation (3.1.2), il est possible de regarder a 
r](x,n) pour x > n > de deux fagons. La premiere fagon est lorsque n est fixe 
et x varie, lequel cas la fonction gi(x) :=: r](x,n) est croissante en x > n. La 
deuxieme fagon, ce qui sera utile, est lorsque x est fixe et n varie, dans ce cas la 
fonction #2 (ft-) :=: i]{ x > n ) es ^ decroissante pour tout n G (0, x). Ainsi pour une 
valeur fixe arbitraire x > 0, il existe n x < x telle que pout tout n avec n < n x , la 
fonction ri~ 1 {x ) n) < 1. Par consequent, 

1 1 

— - < — < 1 pour n <n x 
rj 1 rj 



De meme, pour toute constante k, 

exp ^— j < exp(A;) pour n < n 1 
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De par la remarque 3.1.1 et les equations (3.1.11), (3.1.12), (3.1.22) et (3.1.23), 



: " - - ^, si = ^ et s 2 w 0.010862854400. Si, en plus, n < n x , et, 
en posant les constantes 

5 \ 45 , , / 7 \ 65 



A « s 2 + si + WiSi + exp si + |«j|s 2 + exp ( 



36/ 74 



«2 



+ O.OO83520fexp(O.218423) 



« 0.049784723505 

/ 5 \ 45 
= ^ + ex P(36J74 
« 0.768158487672 



alors 



et 



/ n (z) - cos ^ - ~j 



< 



A B 

- + — 
n r) 



/.(,)= cos (,-I) +0 (I) + (i 



3.1.1. Calculs numeriques des variations totales 

Pour plus de details concernant ce qui suit, se referer principalement a [016]. 
La valeur de Vr(|C|^-Bo(C)) es ^ calculee en trouvant les points stationnaires de 
|C|2.Bo(C) e ^ en utilisant la formule de la variation totale. L'unique point sta- 
tionnaire autre que ( = se situe a ( = C( z i) ~ —1.000459796360 ou Z\ ~ 
1.979495483061. Par consequent, en denotant par 

y(z) = |C(z)|*Bo(CW) 



alors 



Vr(|C| 5 5 (C)) = ^n^W - 2 9( Z = 1) + 2 ^ = *l) + Um 

2— >(J 2— >QO 



1 



+ 2(0.010862854400) + 



12 

0.105059042134 
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Egalement 



lim |C|^o(C) = lim 



lim 



C 



3 5 
' 24 + 24 


- 5 1 


48CU 


3 5 
24 + 24 " 


5 


48C(*)§ 



3 5 

— + 

24 24 



1 

12 

lim |C|^o(C) = lim g{z) 
= 



Pour trouver la valeur numerique de Vr(|C| 2 -£>i(C))) il f &u t l es points stationnaires 
de |C|^i(C). 
En denotant par 

g{z) = \az)\^B l {az)) 
alors ses points stationnaires autres que z — 1 sont 

« 0.138560281581 
2 2 « 1.418538099456 



Par consequent, la variation totale est 

VrOCI^^C)) = lim^)-2^i) + 2^(z = l)-2^ = z 2 )+ lim^) 
« -0.002681327160 - 2(-0.004008186698) 

+ 2(0) - 2(— 0.000639161111) + 
« 0.006613368457 
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avec les valeurs des limites suivantes 

lim |C|^i(C) = lim g(z) 

139 
~ 51840 

lim |C|^i(C) = lim g(z) 



Maintenant, je donne les coefficients A s (() et B s (() qui n'ont pas ete faciles a 
trouver pour s — 1. 

Les coefficients A S (Q et B S (Q pour s = 0, 1 suivent et, dans ce qui suit, t = 

(1 — z 2 )~^, t = —iti avec t\ = (z 2 — 1)~2 et d = —(. 

3* - 5* 3 5 
24(5 "48(2 



B °(() = —777^-77777 C>0 ou 0<^<1 



Si(C) 



Si(C) 



3*1 + 5*? 5 
24Cf 

-(30375* 3 - 369603* 5 + 765765* 7 - 425425* 9 ) 



*o(C) = i -7572 C<0 on z>l 

24Cf 



414720(5 



(405* 2 - 2310* 4 + 1925* 6 ) 

55296C 2 
(1155* - 1925t 3 ) 85085 



110592(2 663552C 5 
(30375*? + 369603*f + 765765*1 + 425425*?) 



C>0 ou < z < 1 



+ 



414720Ci 5 
(405* 2 + 2310* 4 + 1925*j) 
55296Ci 2 



(1155*1 + 1925*?) 85085 

+ ^ + 663552^ C< ° ° U 

110592Ci 2 obd&i)^! 



Ai(C) 
Ai(C) 



81* 2 - 462* 4 + 385* 8 7(3* - 5* 2 



1152 

81* 2 + 462* 4 + 385* 8 
1152 



455 

H52Ci 4608C 3 
7(3*i + 5*?) 455 



C>0 ou 0<^<1 



1152Ci 



+ 



4608Ci 3 



C < ou z>\ 
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Les valeurs de A s (0) et de B s (0) peuvent etre calculees en utilisant les developpe- 
ments de Meissel de J n (n) que nous pouvons retrouver dans [Wa]. Ces valeurs 
se retrouvent egalement dans [016] et, par consequent, 



A 1 (0) 
B (0) 
B 1 (0) 



225 

2l 



140 

(1213.2)1 
~ 2047500 



Remarque 3.1.4 (comparaisons des calculs). Les coefficients mentionnes pre- 
cedemment sont ceux qui ete codes avec Maple pour pouvoir trouver les points 
stationnaires et, par suite, les variations totales. Pour la variation totale de 
Vr(|C| 2-Bo(C))? ^ es t possible de la retrouver avec moins de decimales dans de 
nombreux documents dont [016]. Quant a Vr(|C| 5 -^i(0); e ^ e P eu ^ ^ re retrouvee 
a la page 315 dans [LaWo] oil les auteurs ne donnent pas cependant la fonction 
Bx(C). 

Finalement, les routines en Maple pour evaluer les points stationnaires numeri- 
quement sont dont donnees en annexe. 



3.2. L'ARTICLE DE KUZNETSOV ET FEDOSOV, EXPLICATIONS ET 
CORRECTIONS 

Ayant maintenant une meilleure idee concernant l'expansion uniforme de J n {x), 
alors voici, avec un peu plus details que Particle original [KF], ou il est montre 
que la loi de Weyl pour le disque D de rayon r est dans l'ordre 2/3. 
Sans perte de generalite, le rayon du disque est fixe a r = 1. II est bien connu 
que les valeurs propres du probleme de Dirichlet sur le disque de rayon 1 sont 
les zeros positifs des fonctions de Bessel J n (x) d'ordre n — 0, 1, 2, . . .. Les valeurs 
propres correspondant aux ordres n = 1,2,... ont une multiplicite double et celle 
correspondant a l'ordre n = sont simples. 

Voici une expansion uniforme de la fonction de Bessel de premier type obtenue 
precedemment qui sera utile pour la suite. Soit x > n > 0, J n (x) la fonction de 
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Bessel de premier type et soit 



/„(*) = (l + 6 3 )^(x 2 -n 2 )h n (x) (3.2.1) 

avec 

r)(x, n) = V x 2 — n 2 — n arccos ^— j 

et 5 3 > 0, 5 3 ne depend que de n, \5 3 \ < 0.016456 et Hindoo 5 3 = 0. Alors pour 
une valeur fixe de x, il existe n x < x telle que l'inegalite suivante est verifiee 

f n {x) = cos(^-^)+o(i) + (i) (3.2.2) 

comme il a ete demontre precedemment. 

Remarque 3.2.1. Les nombres de zeros de la fonction f n et J n dans un intervalle 
donne sont les memes, car f n n'est que J n multiplie par une fonction strictement 
positive. 

Remarque 3.2.2. L 'approximation donnee de J n donnee dans Varticle est comme 
suit 

J n {k) = ^-nr^oos^n) - I) +o(I) +o(i)) 

Ainsi, le probleme consistant au denombrement du nombre de zeros positifs des 
fonctions J n (x) n'excedant pas k est transforme en un probleme de denombrement 
de points entiers dans un domaine bien precis d'ou l'utilisation prochaine du 
theoreme de Van der Corput. 

Soit N n (k) le nombre de zeros de la fonction J n (x) dans l'intervalle n < x < k. 
Alors, 

[k] 

N(k) = N Q (k) + J2 N n(k) 

n=l 

ou [k] denote la partie entiere de k. En fixant les nombres z/ = [k^] + \ et 
Vi = [k — &i] + |, la somme N(k) — N (k) est divisee en trois parties comme suit 

[k] v x -\ [fc] 

j2N n (k) = J2 1+ E 1+ E i = Si + s 2 + s 3 . 

Une estimation de S 2 sera donnee en utilisant le theoreme de Van der Corput. 
Quant a et S 3 , leurs estimations utiliseront la monotonicite de N n (k). N (k) 
sera estime a la toute fin. 
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En utilisant (3.2.2) pour n G (uq, Vi) qui est equivalent a aflirmer qu'il existe A 
et B tels que 

A B 

et il s'ensuit que le nombre de zeros de J n (x) pour x G (n, k) est le meme que le 
nombre de zeros de la fonction cos [x — |) dans l'intervalle < x < r](k, n) + — + 
^£^s ■ Par consequent, 



f n (k) - cos (j](k,n) - |j 



< 



1 . 1 A B 
-r]{k,n) + - + - + — — 
7r An r]{k,n) 



Remarque 3.2.3 (Justification du nombre de zeros). Une justification de la for- 
mule pour N n (k) a ete donnee rigourousement par Dominique Rabet. Ses resultats 
paraitront ulterieurement dans ses publications. 

Dans cette partie, une estimation de E 2 est donnee en details. II faut utiliser le 
theoreme de Van der Corput mentionne dans [KF] que voici 
Theoreme 3.2.1 (Van der Corput). Soitles nombres vq—\, — \ etx§—\ entiers 
tel que Vq < V\. Soit la fonction a valeurs reelles f(y) telle que f G C 2 (z/o, V\) et 
telle que 

0<O-<f(u)<T, /»>-, V>l,l2>a- 3 . 

Soit N le nombre de points entiers a Vinterieur du domaine u < v < v\ et 
Xo < x < f(v). Soit A I'aire u domaine. Alors 

N = A + 0(jih) 



Par ce qu'il a ete accompli auparavant, S 2 est donne par le nombre de points 
entiers dans le domaine que voici 

1 ,> v 1 1 A B 
vq<v<v x , - < x < f(v) = -rj(k, v) + - H V— — r. 

2 7r 4 v rj{k,u) 



Soit 5 = max I ^ et a telle que ' 1 ! 

hiy. 



V (k,u) \ < f n Pour z/G (z/ ,^i)- Soit 



11a. 

= -77 + - 5 

ix Am] 

r , \ 1 la. 

h(y) = ~v + - A + — + 6 

7T A TTT) 
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alors fi(v) < f(v) < f 2 (i / )- Par consequent, en denotant N\ et N 2 les nombres 
de points entiers donnes en remplagant / respect ivement par f 1 et / 2 , alors Aq < 
S 2 < N 2 . 

Le theoreme de Van der Corput sera utilise sur fx et / 2 . II faut montrer en details 
que fi et / 2 satisfont aux hypotheses du theoreme. II est suffisant de le faire pour 
/i, le cas de / 2 etant identique. 

Pour montrer que fi satisfait aux hypotheses du theoreme, il faut d'abord evaluer 
les derivees suivantes. 

g = Yv {^/W^-v arccos Q) 



1 



-2v) — arccos ^— j + v 



2 y/k 2 - 1/2 v V/W fl^4 k 

— 1 /z/\ z/ A; 

arccos 



VA; 2 - z/ 2 fc _ z/2 



= - arccos 
<9r? 2 fc 



«9Z/ 2 y/W 



V 



2 



De meme, 



du 7r irr] 2 (k,is) 

= -Vu(k,u)(l+ J 1 

arccosj^ / a 



+ 



et 



7r \ r] 2 (k,u 



d 2 fi _ If -k \ f-f t a \ arccos (f) ^ 2arj u (k,v) 



du 2 n\^k 2 -u 2 J\ V 2 (k,v)J k V V 3 (k,v) 

k 1 / a \ ( arccos (l)) 2a 

7r V/c 2 - v 2 \ V 2 {k,v)J nr] 3 (k,u) 

Comme r](k,v) \ lorsque z/ avec fc fixe et que 77 > 0, alors minr^/c, z/) = 
r](/c, z/x) = 7](/c, k — k^) pour z/ < 1^. L'inegalite arccos(l — a) > y/a pour a £ (0, 1) 
permet done de deduire l'inegalite utile suivante non demontree dans Particle 
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[KF]. 

I k — 

r)(k, k — fcs) = Y k 2 — {k 2 — 2kk^ + ki) — {k — k$) arccos ^ — - — j 



> y 2k^ — k^ — (k — ks)V k~9 



2fcis — A; 18 — k is (Zeis — 1) 

8 l r 18 — 1 1 

V2 



2kw (kw — 2^^) — V^fcis (fcis — 1)- 



> \[2kS 

> k* 

De la meme fagon, 

3 i — 2 2 

r](k : k — k&) > V2/c5 > k$ 
Egalement, l'approximation suivante pour i)(z) = \/l — z 2 — ;zarccos(;z) est utile 

(l-z)l 3 r , s 

!<7= TT < 4^ ^e(o,i) 

V 1 — z — -2 arccos J 4 
Par example, avec < z = r < 1, l'inegalite suivante est verifiee 



\- 2 



^(fc, z/) - -z 2 - ^arccos(z)) 4 A; 

Puisque arccos (y) < arccos (|), alors 

fc~T8 = Vk~^ < arccos(l — k~^) = arccos — arccos ^ 

Par consequent, 

k~w fc _ i8 / a\ <9fi 
< < 1 + — < 



7T 7T \ rf J dv 

De meme, comme < arccos (|) < |, alors 

<9fi 1 /„ o \ 1 a 1 a l,i 

<- 1 + <^ + — T- < ^ + — ttt = ^ + afc" 5 



9z/"2V ^V"2 v 2 {k,k-0) ~ 2 (A;!) 2 2 
II s'ensuit done qu'avec C = - 



< CA;-A < < ^ + 0(A;-5; 



tel qu'ecrit dans Particle [KF]. 



o2 £ 

Pour estimer -g^-, il faut d'abord montrer que 

(arccos(f)) 
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" 'i v") " y \ v — " / 

uniformement pour v G (^0,^1) lorsque — > 00. Pour ce faire, les auteurs separent 
l'intervalle (z/ ,z/i) = (vo,v') U ou u' = k — k^> . Dans l'intervalle (z/ ,z/'), 

comme ^(/c, z/) est toujours decroissante, alors 

(arccos(f)) (arccos(f)) ^ 1 ^m 3 _^ 

7T 2 1 7T 2 1 
. . . < < ; 

_ 4 k ~ 4 _ ^2 

et, dans l'intervalle (i/,!^), 

(arccos(f)) ^ (arccos(^)) 2/ x y ^2/3 1 ^ a 

■■^T(S4-l) H )'-TH( H ) H ) , -'-' t 

Par consequent. 



d 2 A_ 1 / 1 \ g 

Remarque 3.2.4 (Satisfaction des hypotheses du theoreme de Van der Corput). 

Comme 

< Ck~^ < -?A < i + 0(Jfe-5) 
az/ 2 



et que 

du 2 ^ k 

alors f\ satisfait aux hypotheses du theoreme de Van der Corput avec 

a = Ck~^ 

T = l + Oik-*) 

k 
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En denotant ainsi par JVi et A\ respectivement le nombre de points entiers et 
l'aire du domaine delimite par v < v < V\ et x < x < fi(u) et en applicant le 
theoreme de Van der Corput, 

N 1 = A 1 + Ofah) = jH (A(i/) -\)dv + O(fci). 

Comme 5 = max (^), alors 

r bdv < (1/1 - i/o)— < — = 0(A;3). 

./i/ ^0 ^0 

Aussi, puisque r\ v — — arccos (|), que k~Ts < arccos (|) pour v G (fo, ^i), alors 

Ul 1, r lVu , r iKui) dr) r (k ' Ul) dr) 
du — I az/ 



^0 



^ Ju Wu J v (k,u ) J v (k,u () ) - arccos (|)r? 



/, (fc)Vl ) arccos (|)?7 fc-ig \v( k ,^i) . 
Remarque 3.2.5 (correction). Dans V article, une erreur typographique s'est pro- 
duite. En effet, il est ecrit que Ni = 0(k~^$ log A;). 
Ainsi, 

N i = £ (fM-\)du + 0(k-s) 



-J 



^ - -\dv + 0(k*) + 0{k^ log A;) + 0{ki) 

7T 4/ 

f(H)*+«*»> 



II est de meme pour N 2 en montrant de la meme fagon que / 2 satisfait aux 
hypotheses du theoreme de Van der Corput. Par consequent, N-y < £ 2 < N 2 et 

Maintenant pour estimer Si et £ 2 , il faut utiliser la monotonicite de N n (k). De 
fagon generale, 

N n (k) < N n _i{k) pour tout n, k 
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Pour S 3 , comme 

Ik] [k] [k] 



[k-kV]+l 



[fc-fcS]+l 



que 



ri(k,k-k9) 1 A B 

— + 7 + - + - 

7r 4 n T] 



0(ri(k,k-k*)) = 0(k*) 



4 5 

que r)(k, k — fcs) = ^77(1 — et que 



77(1 -ife-i) < (l-(l-jfe-§)) 



fc-t 



4 



r)(k,k — ka) = kfj(l — k 9) 

< A;^ 

alors 

S 3 = 0(£;t£;§) = O(Jfeis) 
Pour Ei, il faut telescoper la somme d'abord comme suit 

"0-5 



Ei = ^iV n (fc) 

n=l 



+ 



+ (^-i(*o-jv§(*)) 
= - 3)^-1 (*) + E - 

n=l 

Maintenant il faut estimer -/V (A;) et N VQ _i(k). Dans le cas de N (k), comme 
J o(^) ~ \/^cos (z - f ), alors iV (A;) = \ + 0(1) et N UQ _i(k) = J + 0(fc§). Par 
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la monotonicite de N n (k), 



N n (k) < N (k) 

N n (k) - N vo _i(k) < N (k) - N Uo _>(k) 

= t + 0(l)-- + 0(kl 

IT IT 

= Oik's) 



et, par consequent, 



2 

n=l 

- (,„-i)(^o(^)) + (, 4)o ( rf) 



Remarque 3.2.6 (correction). Des erreurs typographiques pour I 'expression de 
Ei se sont produites dans V article. En effet, il est ecrit que 



"0-2 "0-2 

£ jv B (fc) = (V - 2) + E - ^0-2 (*)) 

n=l n=l 



devrait s'ecrire 



n=l n=l 
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Ainsi, 

N(k) = N (k) + 2(Ei + S 2 + £ 3 
k 

71 



+ (l) + 2((, D -^ + rt + £'g-i)* + E 



k k C Ul n 2 
- + 2z/ + 2 / ^ - -(ui - v ) + 2E 3 + 0(tf 

7T 7T 7T /,.„ 7T 4 

2z/qA; 



7T 



2 / - -(k - Jfcs - fes) + 2S 3 + O^i 1 

■A* * 2 



-z/ A; + - / 1 ^(fc, - \ + 2S 3 + 0(fc§) 

7T vr ' 2 

2 1 2 Z"^ 1 k 2 f k 2 

-(kv + v O(k-i)) + - Tj(k,u)du-- + - T}(k,v)du + 0(fca) 

r){k,v)dv H — / 7](k,u)du 1 — / r](k, u)du + O(k^) 

~~ vo - - ./ ; . 

2 /"* , . fc 2, 



7T 



/" fc 2 

jf v (k,u)du-- + o(kr 



2(ixk 2 \ k 



ttV 8 / 2 
Ce qui termine l'explication de [KF]. 

Nous rappelons les lecteurs que specialement pour la section 3.2 qui 
s'acheve ici, nous avons utilise k au lieu de A et etudier le probleme 
An + k 2 u = plutot que Au + ku = et cela dans le but de faciliter les 
references a Particle [KF]. Pour la section 3.3 qui suit, nous revenons 
a notre lettre habituelle A et au probleme des valeurs propres exprime 
sous la forme Am + Am = 0. 

3.3. Algorithme pour calculer la fonction de compte 
N(A) 

Cette section est tres courte et comporte une idee que j'ai eue permettant de 
calculer la fonction de compte 



N(X) = 1 
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efficacement sans avoir a calculer toutes les valeurs propres Xj telles que 

A, < A. 

Remarque 3.3.1 (Fagon naive de calculer N(X)). Etant donne I'algorithme de 
la section 2.2 qui determine la suite {Aj}^ =1 pour un entier I donne. La fagon la 
•plus simple de calculer N(X) est de calculer {Aj}'- =1 pour un I suffisamment grand 
et par la suite de compter combien de valeurs propres Xj sont inferieures a A. 
« Buts » des algorithmes 

Etant donne un entier strictement positif /, I'algorithme de la section 2.2 donne 
{Aj}' =1 . Etant donne un nombre reel strictement positif A, I'algorithme de la 
presente section donne l'entier positif N(X). Toutefois de par les remarques pre- 
cedentes, nous pouvons voir qu'il est possible d'utiliser I'algorithme de la section 
2.2 pour evaluer N(X) et cela a un sens en se r appelant les rappels theoriques de 
la section 1.1 concernant N(X). 

Remarque 3.3.2. Les deux algorithmes, celui de la presente section et celui de 
la section 2.2, determinent un ensemble fini de solutions SOL de paires (n,k) 
avec n = 0, 1, 2, . . . et k = 1, 2, . . . ou n et k sont les parametres pour les zeros. 
Les parametres sont I'ordre (n) de la fonction de Bessel et le rang (k) du zero. 
Dans le cas de I'algorithme de la section 2.2, un ensemble fini PA contenant SOL 
est determine a I'aide du theoreme de Courant et SOL est obtenu en ordonnant 
les zeros correspondant aux paires dans PA. Dans le cas de I'algorithme de la 
presente section, seul la frontiere de SOL est evaluee et il n'est pas necessaire de 
trouver tous les zeros qui sont plus petits que \/X pour evaluer N(X). 
L'idee utilise la monotonicite des zeros des fonctions de Bessel J v (x) pour x > 
et v > 0. Si jk(y) denote le k e zero de la fonction de Bessel de J u (x), alors ( voir 
[Wa] ), 

3k{v) < jh{v + v) pour v' > 0. 

Les valeurs propres de A sur le disque D de rayon 1 correspondent aux carrees 
des zeros jk{n) avec k — 1, 2, . . . et n = 0, 1, . . .. En fixant une valeur de A et en 
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denotant par k nt \, la valeur entiere telle que 

jk nA ( n ) < >/A et j kn , x +i{n) > yf\~. 

alors, par la monotonicite des zeros, 

jk n Jn) <j kn+lx (n + l). 

II s'agit done pour chaque ordre n et A choisi a priori d'evaluer k n ^ et ainsi en 
tenant compte de la multiplicity 

N 

N(X) = k Q ,x + 2 k n ,x 

n=l 

OU 

N = max [n ; tel que k n> \ — l}. 

Remarque 3.3.3. L 'ensemble des paires (n,k) solution du probleme denote au- 
paravant par SOL est ainsi done donne par 

SOL={(n,k) , k = l,...,k n , x , n = 0,l,. ..,7V} 

et 

N(X) = 2-card(SOL\{(n,A;) , n = 1, 2, . . . , N, k = 1, 2, . . . , k nA }) 
+ card({(0,A;) , k = 1, 2, . . . , k , x }) 

N 
n=l 

Dans le present algorithme qui calcule N(X), il n'est done pas necessaire d'eva- 
luer tous les zeros associes aux paires dans SOL, mais seulement les valeurs k n> x 
qui determinent la frontiere de SOL. Ce n'est pas le cas pour V algorithme de la 
section 2.2 oil ilfaut tous les evaluer afin de creer {Aj}' =1 et cela s'explique parce 
que les buts poursuivis par les deux algorithmes sont differents quoique mathema- 
tiquement equivalents par les rappels theoriques de la section 1.1. 
L'algorithme en mots 

Les valeurs A^ et A; 0j a peuvent etre trouvees par une methode quelconque. Soit 
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n = 0, 1, . . ., alors k n ^\ — k n +i t \ G {0, 1}. Done pour trouver k n+ i^\ en possedant 
k nt \, il suffit de tester la validite de l'inegalite suivante 

J kn Jn + l)<V\. 

Si l'inegalite est vraie, alors k n+ x t \ = k n \ sinon k n+ i \ = k n \ — l, etc. L'algorithme 
est initialise avec fco,A et l'algorithme procede en reculant et montant jusqu'a N. 
L'algorithme en Maple est annexe a la fin de ce memo ire. 

Remarque 3.3.4 (Precision requise dans l'algorithme). Si nous voulons calculer 
N(X), alors nous voulons trouver tous les zeros j^in) < vX Done par la remarque 
A. 5.1 avec A = 9.00 • 10 08 qui est la valeur maximale de A apparaissant dans le 
tableau 3.1, il faut demander |~log(9.00 • 10 08 )] + 1 = 10 decimales de precision. 
Egalement, le lecteur peut se referer a I 'annexe A. 5 pour plus d 'information. 

3.3.1. Quelques resultats pour les valeurs de la fonction de compte 
N(A) 

Voici quelques resultats pour les valeurs de la fonction de compte. II est interessant 
de constater que pour les valeurs de A = 1 • 10 2 , . . . , 4 - 10 6 dans le tableau qui suit, 
les valeurs de N(X) coincident avec les valeurs retouvees dans la liste des valeurs 
propres calculees par le premier algorithme. Par exemple a l'aide de la liste donnee 
a l'hyperlien precedent, la valeur propre A998978 est la derniere valeur propre plus 
petite que 4.00 • 10 6 et le rang correspond egalement a la valeur apparaissant dans 
le tableau qui suit. 

Voici done quelques valeurs de la fonction de compte. 



57 



A 


N(X) 


A 


N(X) 


A 


N(X) 


1.00E+02 


21 


4.00E+04 


9905 


7.00E+06 


1748690 


2.00E+02 


42 


5.00E+04 


12387 


8.00E+06 


1998600 


3.00E+02 


67 


6.00E+04 


14876 


9.00E+06 


2248481 


4.00E+02 


92 


7.00E+04 


17360 


1.00E+07 


2498404 


5.00E+02 


115 


8.00E+04 


19852 


2.00E+07 


4997783 


6.00E+02 


142 


9.00E+04 


22345 


3.00E+07 


7497255 


7.00E+02 


160 


1.00E+05 


24842 


4.00E+07 


9996841 


8.00E+02 


187 


2.00E+05 


49782 


5.00E+07 


12496477 


9.00E+02 


209 


3.00E+05 


74722 


6.00E+07 


14996075 


1.00E+03 


232 


4.00E+05 


99679 


7.00E+07 


17495807 


2.00E+03 


478 


5.00E+05 


124633 


8.00E+07 


19995483 


3.00E+03 


725 


6.00E+05 


149606 


9.00E+07 


22495362 


4.00E+03 


972 


7.00E+05 


174588 


1.00E+08 


24994959 


5.00E+03 


1214 


8.00E+05 


199546 


2.00E+08 


49992941 


6.00E+03 


1458 


9.00E+05 


224520 


3.00E+08 


74991365 


7.00E+03 


1706 


1.00E+06 


249494 


4.00E+08 


99989990 


8.00E+03 


1952 


2.00E+06 


499298 


5.00E+08 


124988741 


9.00E+03 


2206 


3.00E+06 


749151 


6.00E+08 


149987791 


1.00E+04 


2456 


4.00E+06 


998978 


7.00E+08 


174986827 


2.00E+04 


4925 


5.00E+06 


1248914 


8.00E+08 


199985791 


3.00E+04 


7415 


6.00E+06 


1498757 


9.00E+08 


224984997 



Table 3.1. Valeurs de N(X) pour d ■ 10 p avec d = 1, . . . , 9 et p = 2, . . . , 8. 

3.3.2. Temps de calcul pour les procedures naive et efficace d'eva- 
luer la fonction de compte N(A) 

II est interessant de constater que si notre seul but est d'evaluer N(X), alors il 
est nettement mieux d'utiliser l'algorithme de la presente section et de ne pas 
proceder en utilisant l'algorithme de la section 2.2. Voici a cet effet quelques 
comparaisons des temps de calculs. 



58 



A 


N(X) 


Temps requis (secondes) 


100 


21 


0.024 


1000 


232 


0.118 


10000 


2456 


0.513 


100000 


24842 


1.839 


1000000 


249494 


7.704 



TABLE 3.2. Temps requis pour evaluer N(X) rapidement avec l'al- 
gorithme de la section 3.3 



1 (rang) 


Xi 




Temps requis (secondes) 


21 


98.726272 


103.499453 


0.056 


232 


989.729086 


1001.729979 


0.928 


2456 


9998.868757 


10038.582253 


15.542 


24842 


99997.791736 


100024.359657 


392.066 


249494 


999994.341030 


1000020.429162 


16016.677 



TABLE 3.3. Temps requis pour evaluer N(X) nai'vement avec l'al- 
gorithme de la section 2.2 



Les valeurs de la table 3.2 peuvent etre utlisees pour verifier celles de la table 3.3 
et vice-versa. En d'autres termes, nous pouvons utiliser l'algorithme de la section 
2.2 pour verifier celui de la presente section et vice-versa. 



CONCLUSION 



Finalement, voici une liste resumant les points importants de ce memoire. Les 
points 1,2 et 4 sont des resultats et/ou idees propres a l'auteur de ce memoire. 

(1) La premiere ligne nodale de la deuxieme fonction propre u 2 (r,9) sur un 
secteur S(a) du disque de rayon 1 et d'angle a G (0, 2tt) est donnee par 
1' equation 



r = 



= 



a 



si a < a 



si a > ckq 



ou v = ^ et ou jk{v) est le k e zero de la fonction de Bessel d'ordre v. 
L'angle a est Tangle critique ou la ligne nodale n'est pas definie et il est 
donne implicitement en resolvant l'equation foiy) ~ Ji(2^) = ou v = ^. 
La multiplicite de la 2 e valeur propre vaut 2 lorsque Tangle a = a . 

(2) Soit m G N. Ordonner les valeurs m premieres valeurs propres de A 
sur le disque de rayon 1 ou sur un secteur est equivalent a trouver la 
configuration des ensembles nodaux des m premieres fonctions propres 
qui y sont associees. En utilisant le theoreme de Courant, j'ai cree un 
algorithme permettant de creer la liste des m premieres valeurs propres. 
Une analyse de son temps d'execution en fonction de m pourrait etre fait 
en connaissant une borne sur le temps requis pour evaluer les zeros des 
fonctions de Bessel. 



60 



(3) Soit la fonction de compte N(X) = Ylx <x 1 des valeurs propres du pro- 
bleme An + Xu = avec des conditions de Dirichlet sur 3D. Une explica- 
tion de Particle de [KF] a permis de comprendre pourquoi asymptotique- 
ment 



(4) En utilisant la monotonicite des zeros des fonctions de Bessel ainsi qu'un 
principe de « marche » sur un ensemble de pairs d'entiers bien determine 
par le probleme en soit, j'ai cree un algorithme permettant d'evaluer exac- 
tement et emcacement N(X). Si notre seul but est d'evaluer N(X), le 2 e 
algorithme est beaucoup plus efficace qui si nous procedions nai'vement en 
utilisant le l er algorithme pour trouver le grand index I tel que N(X) < I 
et N(X) > I. 



N(X) 



X 




^ + 0(A*). 



4 



din 
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Annexe A 



LISTE DES PROGRAMMES EN MAPLE ET EN 

MATLAB 

Les programmes donnes ci-dessous ont ete ajuste pour le format de ce memoire. 
Ne faites pas « copier et coller » sans reajuster les lignes des programmes. 

A.l. Programme pour obtenir la suite des m premieres 
valeurs/fonctions propres pour m e N 



# La variable "nbEV" est le nombre de valeurs propres 

# desirees. 
# 

# Le temps requis a la fin de 1' execution est ecrit 

# dans le fichier "ProgState.txt". 
# 

# Thm. de Courant : nb. de comp. conn. <= 0.7*1 ou 

# 1 est le rang. 
# 

# Lancer en arriere-plan: "maple -i Prog2BesFctEV.txt &" 
### 

# The variable "nbEV" indicates the number of wanted 

# eigenvalues . 
# 



A-ii 

# Time required at the of the execution of the program 

# is given in the file "ProgState.txt". 
# 

# Courant Nodal Thm: nb. of connected components <= 0.7*1 

# where 1 is the rank. 
# 

# Send it in background mode : "maple -i Prog2BesFctEV.txt &" 
# 

TimeO:=time() ; 
nbEV: =1000000; 

# nombre de decimal pour les calculs 

# number of decimal for calculations 

# with this accuracy is le-6 for all the zeros. 
Digits :=ceil(log(nbEV)/log(10))+l+6; 

# calcul des grandeurs de tableaux 

# evaluate size of tables 
sizeTAB : =ceil (0 . 7*nbEV) ; 

for nn from 1 to ceil (0 . 5*0 .7*nbEV) 

do sizeTAB : =sizeTAB+ceil (0 . 5*0 . 7*nbEV/nn) 
end; 

TABZero:=Array(l . .sizeTAB); 
TABInd:=Array(l . .sizeTAB) ; 

#cas nn = "ordre" = 



f oo :=1 ; 
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# la mult ipli cite est 1. 

# multiplicity is 1 

for kk from 1 to ceil (0 . 7*nbEV) do 

TABZero [f oo] : =evalf (Bessel JZeros (0 ,kk) ~2) ; 

TABInd[foo] :=0; 

f oo : =f oo+l ; 
end do; 

#STATE 

f idWst : =f open ( " ProgSt ate . txt " , WRITE) ; 
fprintf (f idWst,"\nn = done"); 
f close(f idWst) ; 

#cas nn = "ordre" > 

# la multiplicity vaut 2. 

# multiplicity is 2. 

for nn from 1 to ceil (0. 5*0. 7*nbEV) do 

for kk from 1 to ceil (0 . 5*0 . 7*nbEV/nn) do 
TABZero [f oo] : =evalf (Bessel JZeros (nn,kk) ~2) ; 
#printf ("\n*** %.6f 11 , TABZero [foo] ) ; 
TABInd[foo] :=1; 
foo : =f oo+l ; 
end do; 

##STATE 

#f idWst : =f open ( "ProgState . txt " , WRITE) ; 

#f printf (f idWst , "\nComputation n = °/ d done",nn); 

#f close(f idWst) ; 
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end do; 

# tri sur place 

# sort with respect to a categorical value in right column 

C:=zip((x,y)->[x,y] , TABZero , TABInd) ; 
f : = (x,y)->evalb(x[l]<=y[l]); 
E:=sort(C,f ) ; 
F:=map(x->x[l] ,E) ; 
G:=map(x->x[2] ,E) ; 

#STATE 

f idWst :=f openCProgState.txt" , WRITE) ; 
fprintf (f idWst,"\nSorting done") ; 
f close(f idWst) ; 

TAB_EV : =Array ( 1 . .sizeTAB) ; 

ABC:=1; 

for ii from 1 by 1 while ii<=nbEV do 

#if order=0 then multiplicity = 1; 

if (G[ABC]=0) then TAB_EV[ii] :=F [ABC] ; 

#else order not then multiplicity = 2; 

else TAB_EV[ii] :=F[ABC] ; TAB_EV[ii+l] :=F[ABC] ; ii:=ii+l; 

end if ; 

ABC:=ABC+1; 
end do; 



#STATE 

f idWst : =f open ( "ProgState . txt " , WRITE) ; 
fprintf (fidWst, 11 \nList of EVs created"); 



A-v 

f close (fidWst) ; 

f idW : =f open ( " ListEVs . txt " , WRITE) ; 
for ii from 1 to nbEV do 

fprintf (fidW,'7.d % . 6f \n" , ii ,TAB_EV [ii] ) ; 
end do; 
f close (fidW) ; 

#STATE 

f idWst : =f open ( " ProgSt ate . txt " , WRITE) ; 

fprintf (f idWst,"\nCPU time required = °/,.4f " ,time()-TimeO) ; 
f close(f idWst) ; 

## for computational checks 

#fidW:=f open ("TMP_01.txt" , WRITE) ; 

#for ii from 1 to (foo-1) do 

# fprintf (fidW, 11 °/„d %.6f\n" ,G[ii] ,F[ii] ) ; 

#end do; 

#fclose(fidW) ; 

quit ; 



A. 2. Programme pour evaluer la fonction de compte 



ppBY:=l; 
ppMAX:=10; 
aaBY:=4; 
aaMAX : =9 ; 
abc : =0 ; 

for pp from 1 by ppBY to ppMAX do 
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for aa from 1 by aaBY to aaMAX do 
abc :=abc+l ; 

end do; 
end do; 

SizeSample:=abc; 
SqRLamb:=Array(l . . SizeSample) ; 
TimeTab:=Array(l . .SizeSample); 



Digits:=10; 

# Ajuster boucles en fonction nombre de cas voulus 

# done SizeSample... 
abc : =0 ; 

for pp from 1 by ppBY to ppMAX do 

for aa from 1 by aaBY to aaMAX do 

abc :=abc+l ; 

SqRLamb[abc] :=evalf (sqrt(aa*(10~pp))) ; 
end do; 
end do; 

f idOut :=f open ("ResCount.txt" , WRITE) : 
f close (f idOut) ; 

f idT:=f open ("TimeTrueCount.txt" , WRITE) ; 
f close (f idT) ; 



for aa from 1 to SizeSample do 
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tmpTO :=time() ; 
EigMax : =SqRLamb [aa] ; 

# KO = max(k, zero(ordre 0, index k) <= lambda) 
for KO from 1 to infinity do 

tmp:=evalf (Bessel JZeros (0,K0) ) ; 

if (tmp>EigMax) then kk:=K0-l; break; end; 

end do; 

# Nl = max(n, zero(ordre n, index 1) <= lambda) 
for Nl from to infinity do 

tmp:=evalf (Bessel JZeros (Nl,l)) ; 

if (tmp>EigMax) then sizetab:=Nl; break; end; 

end do; 

Ord:=Array(l. .Nl+1) ; 
Ind:=Array(l. .Nl+1) ; 
0rd[l] :=0; 
Ind[l] :=kk; 

# on monte au ciel 
for nn from 1 to Nl do 

tmp: =evalf (Bessel JZeros (nn,kk) ) ; 

# on recule au Moyen-Age 

for 11 from 1 while tmp>EigMax to infinity do 
kk:=kk-l; 

tmp :=evalf (Bessel JZeros (nn,kk) ) ; 

end; 



0rd[nn+l] :=nn; 
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Ind[nn+1] :=kk; 

if (kk = 0) then break; end; 
end do; 

# Ind[l], c'est pour l'ordre zero done 

# multiplicite est 1 
Tcount:=Ind[l] ; 

for bb from 2 to (Nl+1) do 

# multiplicite est 2 
Tcount : =Tcount+2*Ind [bb] ; 
end do; 

f id0ut:=f open ("ResCount.txt" , APPEND) : 

fprintf (f idOut, 12E %d\n" .round (SqRLamb[aa] ~2) , Tcount) ; 
f close(f idOut) ; 

fidT:=f open ("TimeTrueCount.txt" , APPEND) ; 

fprintf (fidT,"N(°/„.12E) = %d *** Temps requis = °/„.6f 

secondes . \n" , round (SqRLamb [aa] ~2) , Tcount , time () -tmpTO) ; 
f close(f idT) ; 
end do; 



quit ; 
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A. 3. Evaluation des points stationnaires des fonctions 
apparaissant dans le developpement asymptotique 
d'Olver 

Rappelons que 



2 /. t 



2 



< z < 1 



3 



(-Vz 2 — 1 — - arccos f-^ ) 

\2 2 v*// 



2 

z > 1 



alors pour la fonction pour la fonction \((z)\2 B (((z)) ou B est donnee par 

3t - 5t 3 5 
24(2 "48(2 



*o(C) = -^7TT--7^ C>0 ou 0<*<1 



3ti + ht\ 5 
24^ " ~ 48C 2 



fl o(C) = C<0 ou *>1 



et ou 



et ti 



alors voici le code evaluant ses points stationnaires. 



restart ; 
Digits: =20; 

#1 put t below but it IS tl in the text before! 
t:=(z~2 - l)~(-0.5); 

K:=1.5*sqrt(z~2 - 1) - 1 . 5*arccos (1/z) ; 

bO : =(3*t) /(24* (K~ (1/3) ) ) +(5* (t~3) ) / (24* (K~ (1/3) ) ) 



5/(48*(K~(4/3))); 



foo:=diff (((abs(K))~(l/3))*bO,z$l) ; 

ptstat :=f solve (f oo,z=2 .71) ; 

xipt : =evalf (- (subs (z=ptstat ,K) ) ~ (2/3) ) ; 

totval : =(1/12) +2*evalf (subs (z=ptstat, ((abs(K)) ~ (1/3)) *b0)) ; 
fidTerm: =f open (terminal , WRITE) ; 

fprintf (f idTerm, "\ n \n***** Total variation = °/„.20f \n\n" , totval) ; 

f idFile:=fopen("VarTotResultBO.txt" , WRITE) ; 
fprintf (fidFile, "\n\n***** Variation totale = °/„.20f 

\n\n" , totval) ; 

fprintf (fidFile, "\n\n***** Point stationnaire autre que z=l ; 

z = %.20f\n\n" , ptstat) ; 
fprintf (fidFile, "\n\n***** xi (point stationnaire) = 

°/„.20f\n\n\xipt) ; 
fprintf (fidFile, "\n\n***** sqrt ( I xi I )B0(xi) = °/„.20f \n\n" , 

evalf (subs (z=ptstat, ((abs(K))- (1/3)) *b0))) ; 

f close (fidFile) ; 
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Pour la fonction \((z)\2 Bi(((z)) ou B 1 est donnee par 

-(30375* 3 - 369603* 5 + 765765t 7 - 425425t 9 ) 



5i(C) = 



414720(2 
(405t 2 - 2310t 4 + 1925t 6 ) 

55296C 2 
(1155* — 1925t 3 ) 85085 



C>0 ou 0<^<1 



Bi(C) = 



110592^ 663552C 5 
(30375tf + 369603tf + 765765t 7 + 425425*?) 
414720(5 
(405t 2 + 2310^ + 1925tf) 
+ 55296C 2 

(1155ti + 1925tf) 85085 
+ 110592(1 + 663552C C< ° OT 



et ou * et ti sont comme precedemment, alors voici le code evaluant ses points 
stationnaires. 



restart ; 
Digits: =20; 

# 0<z<l ou bien xi>0 

t: = (l-z~2)~(-0.5) ; 

K:=1.5*log((l+sqrt(l-z~2))/z) - 1 . 5*sqrt (l-z~2) ; 

bll:=-(30375*t~3 - 369603*t~5 + 765765*t~7 - 

425425*t~9)/(414720*(K~(l/3))) ; 
bl2:=-(405*t~2 - 2310*t~4 + 1925*t~6) / (55296* (K~ (4/3) )) ; 
bl3 : =- (1155*t-1925*t~3) / (110592* (K~ (7/3) ) ) ; 
bl4:=-85085/(663552*(K~(10/3))) ; 
bl:=bll+bl2+bl3+bl4; 



foo:=diff (((abs(K))~(l/3))*bl,z$l) ; 
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ptstatl : =f solve (foo ,z=0 . 5) ; 

##xipt : =evalf ( (subs (z=ptstat 1 , K) ) ~ (2/3) ) ; 

tmpO : =evalf (limit ( ( (abs (K) ) ~ (1/3) ) *bl , z=0) ) ; 
tmpl : =evalf (subs (z=ptstat 1 , ( (abs (K) ) ~ (1/3) ) *bl) ) ; 

################################################### 

# z>l ou bien xi<0 

tl:=(z~(2)-l)~(-0.5); 

Kl:=1.5*sqrt(z~2 - 1) - 1 .5*arccos(l/z) ; 

blll:=(30375*tl~3 + 369603*tl~5 + 765765*tl~7 + 

425425*tl~9) /(414720* (- (Kl) ~ (1/3) ) ) ; 
bl21:=(405*tl~2 + 2310*tl~4 + 1925*tl~6) / (55296* (Kl~ (4/3) )) ; 
bl31:=(1155*tl+1925*tl~3)/(110592*(-(Kl)~(7/3))) ; 
bl41: =85085/ (663552* (Kl* (10/3))) ; 
bll:=blll+bl21+bl31+bl41; 

fool :=diff (((abs(Kl))* (1/3) )*bll,z$l) ; 

ptstat3 : =f solve (fool ,z=l . 5) ; 

tmp3 : =evalf (subs (z=ptstat3 , ( (abs (Kl) ) ~ (1/3) ) *bll) ) ; 
tmp4 : =evalf (limit ( ( (abs (Kl) ) ~ (1/3) ) *bl 1 , z=inf inity) ) ; 
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totval : =tmp0+tmp4-2*tmpl-2*tmp3 ; 

f idFile:=fopen("VarTotResultBl .txt" , WRITE) ; 

fprintf (f idFile, "\n\n***** Variation totale = °/„.20f 

\n\n" , totval) ; 
fprintf (f idFile, "\n\n***** Point stationnaire no. 1 

autre que z=l ; zl = %.20f \n\n" ,ptstatl) ; 
fprintf (f idFile, "\n\n***** Point stationnaire no. 2 

autre que z=l ; z2 = °/„.20f \n\n" ,ptstat3) ; 
fprintf (f idFile, "\n\n***** sqrt ( I xi I )B1 (xi) pour 

xi(zl) = %.20f\n\n",tmpl) ; 
fprintf (f idFile, "\n\n***** sqrt( |xi I )Bl(xi) pour 

xi(z2) = °/„.20f\n\n",tmp3) ; 
fprintf (f idFile, "\ n \n***** limit ( |xi I )Bl(xi) 

(vers xi(0) = infty) = % . 20f \n\n" ,tmpO) ; 
fprintf (f idFile, "\n\n***** limit ( |xi I )Bl(xi) 

(vers xi(infty) = -infty) = % . 20f \n\n" ,tmp4) ; 
f close (f idFile) ; 



A. 4. Programme en Matlab pour resoudre les problemes 

DE VALEURS PROPRES ET VISUALISER LES ENSEMBLES NO- 
DAUX DES FONCTIONS PROPRES 

Pour de l'information sur ces programmes, vous pouvez taper par exemple 
>>help SpDisque 

dans Matlab dans le repertoire ou vous enregistrez les fichiers pour savoir quoi 
faire. 



function [] = SpDisque (NbEV, NbCourbesNiveau) ; 
% NbEV : nombre de valeurs propres desirees 

°/„ NbCourbesNiveau : quantite des courbes a nouveau pour les figures 
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1 

°/„ Conditions de Dirichlet aux frontieres sont appliquees. 
I 

% DEPENDANCE : DisqueGeo.m 
% 

% Example : SpSecteur (50 , [0 . 25*pi 2*exp(l)] , 1000) ; 
% Conditions aux frontieres 

BdyCond =[1111 1 1 '0' '0' ' 1 ' '0' ; °/„l Dir 
1111 1 1 '0' '0' '1' '0' ; °/„2 Dir 
1111 1 1 '0' '0' '1' '0' ; %3 Dir 
1111 1 1 '0' '0' '1' '0'] ' ; °/A Dir 

°/ Constantes pour la description de 1' equation 
c = 1; a=0; d=l; 

% Efface fichiers information sur mesh et etat du programme 
fidIG = fopenCSpDisqueMeshQual.txt' , 'wt') ; 
f close(f idIG) ; 

fidEt = f open( 'SpDisqueEtat . txt ' , 'wt ' ) ; 
f close (fidEt) ; 

fidEr = f open ( 'SpDisqueErreurs .txt 'wt ') ; 
f close (fidEr) ; 



try 



% Mesh ici 

% [p , e , t] =initmesh ( ' DisqueGeo ' , ' Hmax ' , 2e6 , ' Init ' , ' of f ' , ' Box ' , ' of f ' 

, 'Hgrad' ,1.05, 'Jiggle' , 'minimum' , ' Jigglelter' ,120) ; 
[p, e ,t] =initmesh( 'DisqueGeo' ) ; %%% for testing 
[p, e ,t] =ref inemesh( 'DisqueGeo' ,p, e ,t) ; 
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[p, e ,t] =ref inemesh( 'DisqueGeo' ,p, e ,t) ; 
[p, e ,t] =ref inemesh( 'DisqueGeo' ,p, e ,t) ; 
°/„ [p, e, t] =ref inemesh( 'SecteurGeo ' ,p,e ,t) ; 

°/„ stat sur qualite du mesh 

Mesh_Quality = pdetriq(p.t) ; 

fidIG = fopenCSpDisqueMeshQual.txt' , 'at') ; 

f printf (f idIG, ' \nQualite moyenne des triangles = °/„.6f, 

mean(Mesh_Quality) ) ; 
f printf (fidIG, ' \nEcart-type sur les triangles = °/ .6f, 

std(Mesh_Quality)) ; 
f printf (fidIG, ' \nMinimum de qualite atteint = °/„.6f, 

min(Mesh_Quality) ) ; 
f printf (fidIG, ' \nMaximum de qualite atteint = %.6f\n', 

max(Mesh_Quality)) ; 
f close(f idIG) ; 

% resoudre ici 

MaxEig=1.0; 

while(l) 

fidEt = f open( ' SpDisqueEtat . txt ' , ' wt ' ) ; 

f printf (fidEt ,' \nAppel a pdeeig avec MaxEig = °/„.6f\n\n', 

MaxEig) ; 
f close(f idEt) ; 

[EF_tab,EV_tab]=pdeeig(BdyCond,p,e,t,c,a,d, [0 MaxEig] ) ; 

if length (EV.tab) >= NbEV 
break; 
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else 

MaxEig = 2.00*MaxEig; 
end; 
end; 

°/ fait des beaux dessins ici :-) Sans-Dessin ha ha 

for jj = l:l:NbEV 

fnamestr = char(strcat({'Disque')-,-['EF'}-,-[num2str(jj)}-)) ; 
titlestr = char(strcat({'Disque '},{' \lambda_{'}-, 

{num2str(jj)}, '}=' , {num2str (EV_tab( j j))») ; 
pdeplot(p, [] ,t, 'xydata' ,EF_tab( : , j j) , 'xystyle' , 'flat' , 

'colormap' , 'Jet' , 'colorbar' , 'off') ; 
title(titlestr) ; 

%print ( ' -depsc2 ' , ' -rO ' , fnamestr) ; 
print (' -dj peg' , '-rO' , fnamestr) ; 

end; 

catch 

v = lasterror; 

fidEr = f open( ' SpDisqueErreurs . txt ' , 'at ' ) ; 
f printf (f idEr , '\n°/„s\n' , v. message) ; 
f close (fidEr) ; 

end; 



function [] = SpSecteur (NbEV, angles, NbCourbesNiveau) ; 
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% NbEV : nombre de valeurs propres desirees 

°/ angles : vecteur des angles de chaque secteur 

% NbCourbesNiveau : quantite des courbes a nouveau pour les figures 
% 

% Conditions de Dirichlet aux frontieres sont appliquees. 
% 

% DEPENDANCE : SecteurGeo .m 
I 

% Exemple : SpSecteur(3, [1 .0 1.11 1.112 1.1125 1.11253 1.112531 
1.1125316 1.112531602 1.1125316027 1.11253160273 1.112531602739 
1.112531602739001 1.2], 1000); 

°/ Conditions aux frontieres 

BdyCond =[1111 1 1 '0' '0' ' 1 ' '0' ; °/„l Dir 

1111 1 1 '0' '0' '1' '0' ; V.2 Dir 

1111 1 1 '0' '0' '1' '0'] ' ; °/ 3 Dir 

°/ Constantes pour la description de 1' equation 
c = 1; a=0; d = 1; 

% Efface fichiers information sur mesh et etat du programme 
fidIG = fopenCSpSecteurMeshQual.txt' , 'wt') ; 
f close(f idIG) ; 

fidEt = f open( 'SpSecteurEtat . txt ' , ' wt ' ) ; 
f close (fidEt) ; 

fidEr = f open ( 'SpSecteurErreurs . txt 'wt ') ; 
f close (fidEr) ; 



try 

% Iteration sur les cas differents i.e. les angles 
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for ii = 1 : 1 : length(angles) 

fid = f open( ' AngleGeoSecteur . txt ' , 'wt ' ) ; 
fprintf (f id, '%. 16f\n' , angles (ii) ) ; 
f close(f id) ; 

°/„ Mesh ici 

°/„ [p, e ,t] =initmesh( 'SecteurGeo ' , 'Hmax' , 2e6, ' Init ' , ' of f ' , 'Box' , 

'off ' , 'Hgrad' ,1.05, 'Jiggle' , 'minimum' , ' Jigglelter ' , 120) ; 
[p,e ,t] =initmesh( ' SecteurGeo ') ; %%% for testing 
[p,e,t]=ref inemesh( 'SecteurGeo' ,p,e,t) ; 
[p,e,t]=ref inemesh( 'SecteurGeo' ,p,e,t) ; 
[p,e,t]=ref inemesh( 'SecteurGeo' ,p,e,t) ; 
°/„ [p, e ,t] =ref inemesh( ' SecteurGeo ' ,p,e ,t) ; 

°/„ stat sur qualite du mesh 

Mesh_Quality = pdetriq(p,t) ; 

fidIG = fopenCSpSecteurMeshQual.txt' , 'at') ; 

fprintf (fidIG, '\nSecteur d"angle %• 16f ', angles (ii) ) ; 

fprintf (fidIG, '\nQualite moyenne des triangles = °/„.6f, 

mean(Mesh_Quality) ) ; 
fprintf (fidIG, ' \nEcart-type sur les triangles = %.6f, 

std(Mesh_Quality) ) ; 
fprintf (fidIG, '\nMinimum de qualite atteint = %.6f, 

min(Mesh_Quality) ) ; 
fprintf (fidIG, '\nMaximum de qualite atteint = %.6f\n', 

max(Mesh_Quality) ) ; 
f close(f idIG) ; 



% resoudre ici 
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MaxEig=1.0; 
while(l) 

fidEt = f open( ' SpSecteurEtat . txt ' , ' wt ' ) ; 

fprintf (fidEt , '\nAppel a pdeeig avec MaxEig = °/„.6f\n\n' 

,MaxEig) ; 
f close(f idEt) ; 

[EF_tab,EV_tab]=pdeeig(BdyCond,p,e,t,c,a,d, [0 MaxEig]); 

if length(EV.tab) >= NbEV 

break; 
else 

MaxEig = 2.00*MaxEig; 
end; 
end; 

°/ fait des beaux dessins ici :-) Sans-Dessin ha ha 
for jj = l:l:NbEV 

fnamestr = char (strcat({' Sect '},{num2str(ii)},{'EF'}, 

{num2str(j ; 
titlestr = char(strcat({'Secteur d' 'angle '}, 
{num2str (angles (ii) ,16)}, {' \lambda_{'},{num2str(j j)}, '}=' , 
{num2str(EV_tab(j j))})) ; 

pdeplot (p, [] ,t , 'xydata' ,EF_tab( : , j j ) , 'xystyle ' , 'flat ' , 

' colormap ' , ' Jet ' , ' colorbar ' , ' of f ' ) ; 
title(titlestr) ; 

°/ print ( ' -depsc2 ' , ' -rO ' , fnamestr) ; 
print (' -dj peg' , '-rO' , fnamestr) ; 
end; 
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end; 

catch 

v = lasterror; 

fidEr = f open( ' SpSecteurErreurs . txt ' , ' at ' ) ; 
f print f (fidEr, '\n°/„s\n' , v. message) ; 
f close (fidEr) ; 
end; 

°/„exit ; 



function [x,y] =DisqueGeo(bs ,s) 

°/„ Geometrie du disque de rayon 1 pour resolution PDE 
nbs=4; 

if nargin==0, 

x=nbs; % number of boundary segments 

return 
end; 

d=[ 0.0 0.5*pi 1 0; 11 

0. 5*pi pi 10; °/„2 
pi 1.5*pi 1 0; 13 

1 . 5*pi 2*pi 1 0] ' ; °/A 

bsl=bs(:) ' ; 



if find(bsl<l I bsl>nbs) , 

error( ' Invalid boundary segment.', 
'Non existent boundary segment number. ') 
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end 

if nargin==l, 

x=d( : ,bsl) ; 

return 
end 

x=zeros(size(s)) ; 
y=zeros(size(s)) ; 
[m,n] =size(bs) ; 
if m==l kk n==l, 

bs=bs*ones(size(s)) ; °/„ expand bs 
elseif m~=size(s , 1) I I n~=size(s,2) , 

error( 'PDE:Sector Geometry : SizeBs', 
'bs must be scalar or of same size as s.'); 

end 

if "isempty (s) , 

°/„ boundary segment 1 
ii=f ind(bs==l) ; 
x(ii)=cos(s(ii)) ; 
y(ii)=sin(s(ii)) ; 

°/„ boundary segment 2 
ii=f ind(bs==2) ; 
x(ii)=cos(s(ii)) ; 
y(ii)=sin(s(ii)) ; 

% boundary segment 3 
ii=f ind(bs==3) ; 
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x(ii)=cos(s(ii)) ; 
y(ii)=sin(s(ii)) ; 

% boundary segment 4 
ii=f ind(bs==4) ; 
x(ii)=cos(s(ii)) ; 
y(ii)=sin(s(ii)) ; 

end 



function [x,y]=SecteurGeo(bs,s) 

% Geometrie du secteur de rayon 1 pour PDE solving 
nbs=3 ; 

if nargin==0, 

x=nbs; % number of boundary segments 

return 
end; 

[alpha] = textread( 'AngleGeoSecteur . txt ' ) ; 

if (alpha>=2*pi) 

error( 'BadParam:GeoFile :DesignerCheck' , 

'L' 'angle pour un secteur < 2*pi'); 
end; 



d= [ 0.0 1.0 

. alpha 
1.0 0.0 



1 0; %1 
1 0; °/„2 
1 0]'; V.3 
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bsl=bs(:) ' ; 

if f ind(bsl<l I bsl>nbs) , 

error (' Invalid boundary segment.', 
'Non existent boundary segment number.') 

end 

if nargin==l, 

x=d( : ,bsl) ; 

return 
end 

x=zeros(size(s)) ; 
y=zeros(size(s)) ; 
[m,n] =size(bs) ; 
if m==l && n==l, 

bs=bs*ones(size(s)) ; % expand bs 
elseif m~=size(s,l) I I n~=size(s,2) , 

error( 'PDE:Sector Geometry : SizeBs', 
'bs must be scalar or of same size as s.'); 

end 

if ~isempty(s) , 

% boundary segment 1 
ii=f ind(bs==l) ; 
x(ii)=s(ii) ; 
y(ii)=0.0; 

% boundary segment 2 
ii=f ind(bs==2) ; 
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x(ii)=cos(s(ii)) ; 
y(ii)=sin(s(ii)) ; 

°/„ boundary segment 3 
ii=f ind(bs==3) ; 
x(ii)=s(ii)*cos(alpha) ; 
y(ii)=s(ii)*sin(alpha) ; 

end 



A. 5. Precision des calculs 

Remarque A. 5.1 (Calcul des zeros et precision). Par (1.3.1), la longueur deci- 
male de j u (k) n'excede pas la longueur decimale de v. Par consequent pour calculer 
la partie entiere des zeros jk{v) sans perte de precision, il ne faut que |iog 10 (z^)] 
decimales. Pour etre prudent, nous pouvons exiger |~log 10 (z/)] + 1 decimales de 
precision. Si nous voulons evaluer les zeros a une precision absolue de 10~ p , alors 
il faut que les calculs se fas sent avec |~log 10 (z/)] + 1 +p de decimales. 
La remarque precedente explique pourquoi la ligne de code 
Digits : =ceil (log(nbEV) /log( 10) ) +1+6 ; 

apparait dans le programme de la section 2.2. En effet, le plus petit ordre entier 
non nul n des fonctions de Bessel impliquees dans l'obtention de la suite spectrale 
{\j}"L l est 1. Si, encore, nous denotons par k le rang d'un zero, alors, par le 
theoreme de Courant, 2nk < m comme il a ete explique au chapitre 1 en tenant 
compte de la multiplicite puisque ce sont des valeurs propres de multiplicite 2. Si 
l'ordre est nul, alors k < m et c'est une valeur propre de multiplicite 1 comme il 
a ete explique au chapitre 1. Dans le cas d'un secteur, nk < m. Dans les deux cas 
(le disque ou le secteur), nous avons que l'ordre maximal implique est inferieur a 
m. Puisque m est le nombre de valeurs propres desirees c'est-a-dire m :=: nbEV 
dans le programme, alors nous voyons d'ou sort la formule. Le choix d'imposer 
une precision de 10~ 6 s'explique par le fait que nous voulons eviter de trouver 
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des zeros ayant la meme representation jusqu'a 1CT 3 comme il est possible d'en 
trouver dans la liste de 30 MB des valeurs propres Xj pour j — 1, . . . , 10 6 . 
Dans le cas du programme de la section 3.3, il faut evaluer tous les zeros jk(n) < 
Puisque la valeur maximale de A pour laquelle nous voulons evaluer N(X) 
est 9.00 • 10 s , alors le nombre de decimales requises [log 10 (9.00 • 10 8 )] + 1 + = 10 
d'ou la ligne de code 
Digits :=10; 

dans le programme de la section 3.3 evaluant N(X). 

A. 6. Methodes utilisees par le logiciel Maple pour eva- 
luer LES ZEROS 

Ci-dessous, nous trouvons une bonne partie de l'information necessaire pour savoir 
comment Maple evalue les zeros des fonctions de Bessel. C'est un courriel qui m'a 
ete envoye par l'equipe technique de Maplesoft. 
Comments extracted from the source file: 

Rayleigh's formula for the first zero of Bessel J function can be 
obtained in reference: 

G.N. Watson - A Treatise on the Theory of Bessel Functions, 
2nd ed., CUP, 1966, page 502, 

Olver's asymptotic expansion for large order can be obtained in 
references : 

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, 
Dover, NY (1972) 

Royal Society Mathmatical Tables, vol. 7, part III. Ed. F.W.J. 
Olver, CUP, 1960 
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Olver, F.W.J. , Philosophical Transactions of the Royal Society, A, 
vol 247, pag. 307-368 (1954) 

McMahon's asymptotic expansion for large n can be obtained in 
references : 

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, 
Dover, NY (1972) 

G.N. Watson - A Treatise on the Theory of Bessel Functions, 
2nd ed. , CUP, 1966 

To find the N_th ZERO OF BESSEL FUNCTION of order m, 
the procedure uses the following results proved by 
Schafheitlin see 

G.N. Watson - A Treatise on the Theory of Bessel Functions, 
2nd ed., CUP, 1966: 

For m<=l/2 the n_th zero of BesselJ [m] (x) is in the interval: 
[(n+m/2-l/4)*Pi .. (n+m/4-1/8) *Pi] 

For l/2<m<5/2 the n_th zero of BesselJ [m] (x) is in the interval: 
[(n+m/4-l/8)*Pi .. (n+m/2-1/4) *Pi] 

For m>l/2, the zeros of Bessel J [m] (x) which exceed 
(2*m+l)*(2*m+3)/Pi are in the intervals: [(p+m/2+1/2) *Pi .. 
(p+m/2+3/4) *Pi] where p are consecutive positive integers. The 
relation of p with n has been taken as p=n-l. 
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The zeros less than (2*m+l) * (2*m+3) /Pi will be found by a recursive 
algorithm. 

We use the following assymptotic aproximations for the first and 
second zeros [see Watson]: 

first zero: m + 1 .85*m~(l/3) ) + 0(m~(-l/3)) 
second zero: m + 3 . 24*m~ (1/3) ) + 0(m~(-l/3) ) 

To find the n-th zero j(m,n) we use the interval 

[j(m,n-l) .. j(m,n-l) + incr] 
where the increment is 

incr = j(m,n-l) - j(m,n-2) 

The other part of the solution can be found by setting a Maple 
Interface variable to display the initial routine that is used. 

interf ace(verboseproc=3) ; eval (BesselJZeros) ; 

I hope this helps answer your questions. 

Cheers , 

Ian Maplesoft Technical Support www.maplesoft.com/support/ 
Original Message 

From: Maplesoft Customer Service Sent: Monday, November 05, 2007 
4:15 PM To: Claude Gravel; Maplesoft Technical Support Subject: RE: 
Customer Service Request: (Web) Other 
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Comme nous pouvons le constater, Maple utilise les developpements asympto- 
tiques d'Olver pour representer les fonctions de Bessel pour des ordres grands. 
Les developpements asymptotiques d'Olver des zeros pour des ordres grands sont 
egalement utilises pour evaluer les zeros. 



